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Introduction 



1 Regulations on The Olympiad 

Regulations on holding The Olympiad contest for students on Mathematical 
and Theoretical Physics were developed in April 2010 [see Special Issue No. 2]. 
They were signed by the three parties: Samara State University (hereinafter referred 
to as SamGU), Steklov Mathematical Institute (SMI RAS), and Moscow Institute 
of Physics and Technology (MIPT). The text of the regulations is given below. 

2 Carrying out The Olympiad 

On May 21-24th, 2010, Ail-Russian Student Training Olympiad in Mathematical 
and Theoretical Physics "Mathematical Physics" with International Participation 
has been held. It was the second in the series of Olympiads. It is planned that in 
future such Olympiads will take place annually. 

The organizers of the series of Olympiads on Mathematical & Theoretical Physics 
"Mathematical Physics" are: 

Aleksander Anatolyevitch Andreyev (staff member of the Scientific Research 
Laboratory of Mathematical Physics of SamGU), 

Georgiy Sergeyevitch Beloglazov (The University of Dodoma- UDOM, Tanzania; 
Perm State Pharmaceutical Academy), 

Boris Vasilyevitch Danilyuk (staff member of the Scientific Research Laboratory 
of Mathematical Physics of SamGU), 

Mikhail Vyacheslavovitch Dolgopolov (Head of the Scientific Research Laborato- 
ry of Mathematical Physics of SamGU), 

Vitaliy Petrovitch Garkin (Vice-rector for Academic Affairs, Chairman of the 
Local Organizing Committee, SamGU), 

Mikhail Gennadievich Ivanov (Associate Professor, MIPT), 

Yuri Nikolayevitch Radayev (staff member of the Scientific Research Laboratory 
of Mathematical Physics of SamGU), 

Irina Nikolayevna Rodionova (staff member of the Scientific Research Laboratory 
of Mathematical Physics of SamGU), 

Yuri Aleksandrovitch Samarskiy (Deputy Vice Chancellor on Education, MIPT), 

Irina Semyonovna Tsirova (docent, SamGU), 
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Igor Vasilyevitch Volovich (scientific leader of the Scientific Research Laborato- 
ry of Mathematical Physics of SamGU, head of the department of Mathematical 
Physics of MIAN), 

Aleksander Petrovitch Zubarev (staff member of the Scientific Research Labora- 
tory of Mathematical Physics of SamGU). 

The Olympiad has been held as a team competition. Number of participants 
of each team - from 3 to 10 students of 2nd to 6th courses (years) of higher educa- 
tional establishments of Russia, CIS, and other countries. It was allowed that more 
than one team participates on behalf of any organization. Order of the Olympiad: 

The participants have been offered to solve 14 problems. Time to start solving 
problems of the contest was 11:00 pm Moscow time on May 20th, 2010. The state- 
ments of the contest tasks are published in *.pdf format at the webpage of the 
Olympiad www.labmathphys.samsu.ru/eng/content/view/29/36/ of the website of 
the Scientific Research Laboratory of Mathematical Physics of SamGU 

www . labmat hphy s . samsu . ru / eng 

and have been sent to the registered participants of The Olympiad. 

The deadline to send the scanned (or photographed) solutions to the E-mail 
address of the Mathematical Physics Laboratory: slmp@ssu.samara.ru was 11pm 
Moscow time on May 24th, 2010. All participants of The Olympiad who had sent 
their solutions by E-mail, have received confirmation that their solutions had been 
accepted. 

It was allowed that the participants solve any problems from the number of the 
proposed ones which they find affordable for the own level of knowledge digestion in 
different units of mathematics and physics thus participating in the topical scoring 
nomination (for purpose of this scoring nomination, the problems are aggregated 
into groups 1 to 3 problems in each). 

In the application letter, the name of organization hosting the team should be 
stated together with the surname, name, for each participant of the team, Depart- 
ment (speciality), course/year; contact E-mail address. 

The Nominations of The Olympiad: 

1) The overall team scoring based on the three best team participants perfor- 
mance (3 prize-winning team places). In the present Olympiad, it is possible to 
submit only one solution on behalf of a team; it is advised to mention the author(s) 
of every solution or solution method [stating also the year(s)/course(s) of studying] 
at the end of each solution (or method of solution). The winner is the team which 
participants have solved correctly maximum number of different problems. Any par- 
ticipant of a team has the right to send a solution separately. Within the team 
scoring, the correct solutions will be considered and accounted. The maximum pos- 
sible number of points in a team scoring is 14 (because the total number of problems 
offered is 14). 

2) It is possible for a student to participate in the overall personal contest (with- 
in the framework of the Olympiad by correspondence) ON CONDITION OF THE 
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PRESENCE OF AN INDIVIDUAL APPLICATION (REQUEST) from a partici- 
pant of The Olympiad (3 prize-winning places). 

3) Overall team topic scoring (1-3 prize- winning places on each subject). 

4) Separate team scoring among each of the years (second through sixth courses). 

5) Best team among the technical specialities of the institutes of higher education. 

6) Other nominations. Separate nomination is supported by the Center on Ad- 
vanced Training and Professional Development at Samara State University. 

The winners of The Olympiad held by correspondence participated in the day 
competition Olympiad held in Samara in September - 2010 (at the same time with 
the Second International Conference and School on Mathematical Physics and its 
Applications). For the above said winners, their travel and/or accommodation ex- 
penses were reimbursed. 

3 Contents of the problems for the Olympiad contest 

The topic range of our 'Olympiad' is related to mathematical methods in describ- 
ing physical phenomena based on the following units of mathematics and theoretical 
physics: 

theory of differential, integral equations, and boundary-value problems; 
theory of generalized functions, integral transform, theory of functions of complex 
variable; 

functional analysis, operational calculus, spectral analysis; 
probability theory, theory of random processes; 
differential geometry and topology; 

theoretical mechanics, electrodynamics, relativity theory, quantum mechanics, 
and gravitation theory. 

New scientific methodological approach to composing the statements of the prob- 
lems for The Olympiad was first introduced in the sense that about a half of the 
problems offered to the participants for the solution supposed that certain stage of 
research (taken from original modern academic research in mathematical physics 
and its applications) is involved. On the basis of the above mentioned approach, the 
recommendations on composing statements of the problems for The Olympiad are 
developed. 

In the present issue, we quote the statements of problems offered to the par- 
ticipants of All-Russia Students Training Olympiad in Mathematical and Theoreti- 
cal Physics "Mathematical Physics" with International Participation (held on May 
21-24th, 2010). 

4 Results and resume of The Olympiad 

In The Olympiad, the teams from the following institutes of higher education 
and other organizations have participated: 
Belarusian State university, 
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Moscow Institute of Physics and Technology (State University), 
National University of Singapore, 

Department of Theor. Phys. named after I.E.Tamm of FIAN (the Institute of 
Physics of Academy of Sciences of Russia) - postgraduate, 

Samara State University of Architecture and Construction (two teams), 

Samara State Aerospace University (SGAU), 

Samara State University, 

The Federal University of Siberia, 

Ulyanovsk State Pedagogical University, U1GPU (two teams), 
The University of Dodoma (UDOM, Tanzania), 
Yaroslavl State University (the team of the Physics Department). 
The jury has positively assessed the works by the following participants of the 
teams: 

Belarusian State University: Alexey Bobrick. 

Moscow Institute of Physics and Technology (State University): Kostjukevich 
Yury, and the fourth course team: Nikolai Fedotov, Anton Fetisov, Mikhail Shalagi- 
nov, Aleksander Shtyk. 

Department of Theor. Phys. named after I.E.Tamm of FIAN (the Institute of 
Physics of Academy of Sciences of Russia): Andrey Borisov. 

Samara State University of Architecture and Construction, SGASU (two teams 
of the students of the 5th year): leader - Sergey Zinakov. 

Samara State Aerospace University: Mikhail Malyshev, Yekaterina Pudikova. 

Samara State University: team of theoretical physicists - Tatiana Volkova, Matvei 
Mashchenko, Maksim Nefedov, Yelena Petrova. 

The Federal University of Siberia: Artyom Ryasik, Polina Syomina, Anton Sheykin. 

Ulyanovsk State Pedagogical University (two teams): 3rd year - Yuri Antonov, 
Aleksandra Volkova, Oksana Rodionova; 

5th year: Maria Vasina, Artyom Ovchinnikov. Aleksander Chaadayev, Aleksander 
Ernezaks. 

The diploma of Laureates or diploma of the winners in nominations have been 
sent to all above mentioned participants. All participants of The Olympiad have 
been invited to attend the School-2010 on Applied Mathematical Physics (PMF) 
from July 1st till July 14th, and the scientific Conference together with another 
School & Olympiad (August 29th - September 9th, 2010). 

The Winners of The Olympiad in the nominations: 

The Overall Team Score: 

1st Place, 5 problems solved correctly (means, 8 and more points per a problem, 
maximum 10 points per a problem), the total score is 102 points, - the team of the 
4th year of MIPT. The winners are granted the prize - traveling costs be paid for 
them to participate in the scientific Conference together with School & Olympiad 
(August 29th - September 7th, 2010). 

2nd Place - The Federal University of Siberia 

3rd Place - Samara State University 



8 



Modern Problems of Mathematical Physics. Special Issue N» 3 



The Total Personal Score: 

1st Place, 4 problems solved correctly, score is 126 points, - Alexey Bobrick, 
Theoretical Physics magistracy at the Department of Physics of Belarusian State 
University. The winner is granted the prize - either traveling or accommodation costs 
be paid for him to participate in the scientific Conference together with School & 
Olympiad (August 29th - September 7th, 2010). 

2nd Place, 3 problems solved correctly, score is 74 points, - Yuri Kostyukevitch, 
the student of the 5th year of the Department of Molecular and Biological Physics, 
group No. 541, MIPT. Recommended for the Magistracy or (post)graduate school 
of MIPT. 

3rd Place, 2 problems solved correctly, score is 70 points, - Anton Sheykin, the 
student of the 4th year at Engineering Physical Department of IIFiRE (Physics and 
Radioelectronics) of Siberian Federal University. Recommended for the Magistracy 
or (post)graduate school of MIPT. 

The best (complete and original) solutions of separate problems: by Alexey Bo- 
brik, Polina Syomina, Mikhail Shalaginov, Anton Sheykin. 

1st Team Place among the 3rd year students - SamGU; 

2nd Team Place among the 3rd year students - SGAU. 

1st Place in Personal contest among the 3rd year students - Maksim Nefedov; 
2nd Place in Personal contest among the students of the 3rd year students - 
Mikhail Malyshev. 

1st Place in Personal contest among the 4th year students - Anton Sheykin; 

2-3 Places in Personal contest among the students of the 4th year - Nikolai 
Fedotov, Anton Fetisov, and Mikhail Shalaginov. 

Among the teams of Pedagogical, Engineering & Technical institutes of higher 
education: 

1st Place - U1GPU, 3 year; 

2nd Place - U1GPU, 5 year; 

3rd Place - SGASU. 

All winners and prize winners of The Olympiad are granted with the free of 
charge accommodation at the PMF School-2010. 

The information on the Olympiad, formulation of the Problems-2010 statements, 
answers and solutions of the tasks-2010 are presented in this document. 

www . lab mat hphys . samsu . ru / eng simp @ ssu . samara, ru 

Organizers of a series of the Mathematical Physics Olympiads: Alexander Andreev, 
George Beloglazov, Boris Danilyuk, Mikhail Dolgopolov, Vitaliy Garkin, Mikhail Ivanov, 
Yury Radaev, Irina Rodionova, Yury Samarsky, Irina Tsirova, Igor Volovich, 

Alexander Zubarev 



Problems and Solutions 




1. Virial for anharmonic oscillations 

For a particle moving along the x axis with Hamiltonian 



where A is a positive constant, m is the mass of the particle, p is the momentum of 
the particle, n — 1, 2, 3, ... , obtain the relationship between the average values 
of kinetic (K) and potential energy (U) using two methods: 

(a) directly from the virial theorem (see explanation below); 

(b) from the condition 



which is true due to the fact that the motion of the particle is finite. 

Instruction : when a particle moves in a potential field, its Hamiltonian H and 




acting force F are defined by 



H = K + U, 



F = — grad U. 



Explanation to Problem 1. In classical mechanics time average values of kinetic 
and potential energies of the systems performing finite motion are in rather simple 
relationship. 
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The average value for a physical quantity G for a sufficiently large time interval 
r is defined in a standard way: 



T 

(g) = i y g dt. 



If (if) is the average (for a rather long time interval) kinetic energy of the system 
of point particles (radius- vectors of the particles given as fi) subjected to forces Fi, 
then the following relation takes place: 

The right hand side of equation (♦) is called Clausius virial, and the equation itself 
expresses the so called the virial Theorem. The proof of the theorem is given, for 
example, in [T]. 

SOLUTION 

(a) According to the given statement the particle is moving in the field of a 
potential force and possesses potential energy U(x) = Ax 2 ™. The force equals 

Fx = JJL = -2\nx 2n -\ 
ax 

Substitute it into the equation (♦) which expresses the virial theorem: 

(K) = --(F ■ r) = --(F x x) = ]^{2\nx 2n ~ l x) = n(Xx 2n ) = n(U). 

(b) According to the given conditions, 



2n-V 



I d \ I dx dp \ / P \ I p 2 

0= \dt {xp) / = \lI P + X Tt/ = \m P + XF */ = \m + Xi ~ X2nX 

2 

= ( — )- 2n(\x 2n ) = 2(K) - 2n(U), 
\m/ 

that is why 

(K)=n(U). 

2. Method of successive approximations 

Solve the integral Volterra equation of 2nd kind 

v(x) = T ^^ ) Jv(t)dt + f(x), (1) 
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where x G [0,h], f(x) is a given (known) continuous on [0, h] function, a(x) G 
C 1 [0, h] (continuously differentiable function), and a(x) 7^ 1, a'{x) is the derivative. 
Perform your solution check. 



Vito Volterra (3 May 1860 - 11 October 
1940) was an Italian mathematician and 
physicist, known for his contributions to 
mathematical biology and integral equa- 
tions. 



SOLUTION 

Using method of successive approximations, we find the solution of the equa- 

a'(x) 




tion ([T]) through kernel resolvent K(x, t) 



1 — a(x) 
cp{x) = I R(x,t)f(t)dt + f(x), 



(2) 



R(x,t) = J2 K n(x,t), 



n=l 



ATi(x,t) = K(x,t), 

X 

K n (x,t) = j K 1 (x,s)K n - 1 (s,t)ds. 
t 

Using formula (jSJ) we find the repeated kernel K 2 (x,t), K 3 (x,t) 

X 

a'(x) a'(s) 



K 2 (x,t) 



1 — a(x) 1 — a(s) 



ds 



n'(r) 

— { -J-[Hl-a(t))-Hl-a(x))], 



(3) 
(4) 
(5) 



(6) 



K 3 {x,t) 



a (si 



a'(x) 

1 — a(x) J 1 — a(s 
t 



[ln(l -a{t)) -ln(l - a{s))]ds 



(7) 



12 



Modern Problems of Mathematical Physics. Special Issue N» 3 



a'(x) [ln 2 (l — a(t)) , .„ .„ ... ln 2 (l — a(x)) 
— — - ln(l - a(t)) ln(l - a(x)) + — — 



1 — a(x) 



a'{x) 



2!(1 -a{x)) 
Similarly, using the formula (jSJ) 

a'(x) 



[ln(l -a(t)) -ln(l - a(x))] 2 . 



# 4 (M) 



3!(1 -a(x)) 
we come to the conclusion that 

a'(x) 



[ln(l -a(t)) -ln(l - a(x))f 



K n (x,t) 



(n-l)!(l-a(x)) 



ln r 



! / 1 - a(t) 
1 — a(x) 



(8) 



(9) 



Expression (jUJ) should be substituted into the formula 

' l-q(t) 

1— 



J R(x,t) = ^ T 



ln n-l f l-a(t) 

a (x) 



n=l 



a 



(x) (n-1)! 



(10) 



Performing the transformation n — 1 = m in (fTDjl and recalling the expansion 



z — ' m! 



m=0 



we obtain 



R(x,t) 



a'(x)(l - a(t)) 



(l-a(x)) 2 

Substituting (1121) into formula ([21 , we obtain 



(11) 



(12) 



<^(x) = /(x) + 



a'(x) 



/(t)[l-a(t)]dt. 



(13) 



[l-a{x)f 
o 

Checking. Let us show that the function (I13p is the solution of eq. ([I]). Designate 

X 

a'(x) 



J(x) = <f(x) 



1 — a(x) 



ip(t)dt, 



(14) 



and substitute the function (TT31 into the right side of equation (I14p . As a result, we 
shall obtain 



J(x) = /(x) + 



a'(x) 
[1 - a(*)] s 



/(t)[l-a(t)]dt- 



(15) 



Mathematical Physics: Problems and Solutions 



13 



a'(x) f a'(x) f f a'(t)(l — a(s))f(s)ds 



1 - a{x) J 1 - a(x) J J [1 - a(t)]' 

o oo 



In the last term of formula (115]) let us change the integration order and calculate 
the inner integral: 



a'(t)dt 1 1 



[l-a(t)] 2 l-a(x) l-a(s)' 

s 

The result should be substituted into the formula (flBl: 



(16) 



J(x) = f(x) + p 5^i)]2 / /(<)[! - «(<)]*- (17) 



f f(t)dt- H a ' {x )^ f(l-a(s))f(s)ds + -^^ 



1 — a(x) y (1 — a(x)) 2 J 1 — a(x 

o o o 

Our checking has shown that the function (TIB"]) is the correct solution of the 
equation ([TJ. 

3. Evaluation for ultrametric diffusion 

When solving equations of the ultrametric diffusion type (that have a relation to 
the description of conformational dynamics of complicated systems such as biomacro- 
molecules) the results can often be presented in the form of series of exponents. Two 
of such series are represented below: 



R(t) = a- n e- b ~ n \ S(t) = ^ a ~ n ^ nt - 

n=0 n=l 

Here t is time, R(t) and S(t) are probabilities that a system is in some definite 
groups of states, k is some integer number, a > 1, b > 1 are some parameters. 

Study the asymptotic behavior of functions R(t) and S(t) at t — > oo and evalu- 
ate, if possible, their asymptotics using elementary functions depending on t. 

SOLUTION 

Let us explore S(t) and R(t) = 5(^1^=0 + e~*. Note that the function — a~ x 

decreases, while the function e~ b * increases with the growth of x. Then in the 
interval x — 1 < n < x the inequality takes place 

io-V»" (,_1)< < - h a- n e- b ~ nt < —}—-a-^e- b - Xi 
x k n k (x — \) k 
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takes place. Integrating it with respect to x from n to n + 1 gives (for n > 1): 

-71+1 1 1 /-n+1 j 



in % n In 

Now, by summing over n from 2 to oo, we obtain: 

S m - m (t) ^S(t)^ S max (t) 

where S(t) = ~, 



(x-r 



a- x e- bn dx. 



h a- n e- b ~ n \ 

2 



roc -i roc -i 

S min (t) = a" 1 / —a-^e~ h ^ lh dx, S max (t) = a / —a^e^cfe. 

J2 X J2 \ X ~ J-J 

By switching to new variables, we have: 

S mm (t) = a" 1 (hi ft)*" 1 (lnt)-V£f jf* (l - ^Q^' y^e^dy, 



5 max (t) = aflat)*- 1 (In j* '* (l - ^y^fef-Vdy. 

Let us designate the function 

7 W (*, * , «, /?) - (l - ^) ^ y^dy, 

considering a (3 < 1 for convergence. 

Note that the limit for this function is the Gamma-function (see the proof in j2]): 

rOO 

lim 7 (fc) (z, t, a, (3) = / y z ~ l e~ y dy = T{z). 
Then at t 3> 1 it is possible to write 

S mm (t) = a- 1 (Inb)"- 1 (IntyT^r (j^) (1 + o(t)) , 

S max (t) = a (lnb) k ~ l (lnf)~* r&r (1 + o(t)) . 

Notation o(t) means that in the limit at t — > oo the value of o(t) tends to zero. 
Since 

S(t) = S(t) + a- 1 e- b ~ lt , 

In a 

the final asymptotic evaluation is of the form (In t) ^ t In b : 



In a 
mT 



(lni)- fc r^f(l + o(t)) < 



Mathematical Physics: Problems and Solutions 



15 



<a(ln6) fc - 1 r(^)(lnt)- fc rfef(l + o(t)) 



Because 

R(t) = S{t)\ k=0 + e~ t + a~ 1 e~ h ' 1 \ 

we have also 

a-^ln&^r r^s (1 + o(t)) < R(t) < aQn&^r 



In a 
\nb 



In a 



(i + (t)) 



4. Double effort 

Solve the Volterra integral equation 

<p(x) = x + (s — x)ip(s)ds. (1) 








SOLUTION 

The considered equation is in fact Volterra integral equation of 2 nd kind with 
continuous kernel. According to the theory of linear integral equations, it has a 
unique solution. To find its solution, we reduce it to Cauchy problem for ordinary 
differential equation. Q 

Let us assume that <p(x) is the solution of the equation (CQ). Double differentiation 
of the identity equation §1§ by x gives: 



<pf(x) = l + / (-tp(s))ds, (2) 
Jo 

<f(x) = -<p{x). (3) 

<p(x) = C\ sinx + C*2 cos x. (4) 
Assuming x = results in: 

V(0) = 0. (5) 



1 Special Issue No. 4 contains also another methods of solution (see [2]). 
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At x = we have from eq. ([2]): 

</(0) = 1. (6) 

Hence C\ — 1, C 2 = 0, and 

= sin re. (7) 

It is possible to perform a check: 

rx 

J(x)=x + / (s — x)smsds= (8) 
Jo 

/•a; 

= x — (s — x) cos s\q + / cos s ds, (9) 

= x — x + sin(s)|g = sinx. 

It is possible to propose other ways to solve this equation, such as with the help 
of Laplace transformation, or by building a resolvent of the kernel by successive 
approximations method. However, the technique developed above is the simplest. 



5. Random walk 

A particle performs random walk on one-dimensional lattice situated on the 
OX axis, the nodes of the lattice have the coordinates m = 0, ±1,±2, .... At the 
initial time moment t = the particle is at the origin of the coordinates. At random 
time moments ti, t^ t%,. . . the particle performs the jumps into adjacent lattice nodes 
with the probabilities of a jump leftwards and rightwards equal to — , the probability 

of remaining still being — 1 — a. The time intervals between the jumps t i+ i — ti, 
i = 0, 1,2, ... are independent random quantities which have the same exponential 

distribution = — exp(— t/r) with expectation r. Find: 

r 

(a) dispersion of the location of the particle as time function t; 

(b) probability that the particle is in m-th node at time moment t. 

SOLUTION 

1. Let us find the probability p(t,n) that the particle during time interval (0, i] 
would perform exactly n jumps, taking into account that this probability is the 
distribution function of the Poisson process. Let tx,t2,-.. be time instants of the 
jumps, so that 

t Q = < h < t 2 < < t„_l < t„ < t < t n+l . (1) 

Probability p(t, n) that the particle during the time interval (0, t] would perform 
exactly n jumps, can be represented as 



p(t,n) = M[I {t n < t < t n+1 )) , 



(2) 



Mathematical Physics: Problems and Solutions 



17 



where M [...] is expectation and 
I(t n <t< t n+1 ) = 
Let us perform Laplace transformation of the function p(t,n): 



1, if t n < t < t n+ i, 

0, if t n > t or t > t n+ \. 



p{s,n) 



/ dte~ st p(t, n) = M / dte~ st I (t n < t < t n+1 ] 
Jo Uo 



(3) 



M 



g S^n g -^n + 1 



(4) 



Due to the fact that t n = r» is the sum of independent random variables, then 



i=i 



exp — s ) 

i=l 



]M[exp(-sr i )] 



(5) 



i=l 



n roo 

[ / rfr J exp(- S r i )$(r J ) = <l n (s), 



where $(s) = / rfrexp (— sr) $(t) is Laplace image of $(r). From this, 



p(s, n) = M 



1 - $fs] 



(6) 



Due to the fact that = — exp(— t/r), it follows that 

r 



and 



p(s,n) 



1/r 



s + l/r' 

(1/t) w 
(s + l/r) n+1 



(7) 



(8) 



Transforming from Laplace image to the original, we obtain the distribution function 
for the Poisson process 



p(t,n) 



t n e -t/r 



T n n\ 



(9) 



2. Let us find the dispersion D(t) of the location of the particle as a function 
of time t. Location of the particle after n jumps is defined by the random variable 
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X n (t) = £i + £ 2 + ••• + ^n, where % — 1, are independent random variables, 
possessing the values ±1 with the probability — and with the probability j3 = 

= 1 - a. Due to the fact that M [(&) 2 ] = a, M(&) = and M = for i ^ j, 
the dispersion equals to 

oo oo 

D(t) = J>(*,n)M [(Xn(t)) 2 } = £>(t,n)cm. (10) 

n=0 n=0 

Substituting into this formula the expression for p(t,n), we obtain 

D (t) = Y i^—an = ae-^- V ± = a-. (11) 

v ' ^ r n n\ t ^ T n ~ x (n-l)\ t 

n=l n=l v ' 

3. Let us find the probability that the particle is located at node m at time 
moment t. Let h n (m) designate the probability that the particle is located at the 
point with coordinate m after n jumps (transitions). Then the probability f(m,t) 
that the particle is at location m at time moment t will be given by formula 

oo 

f{m,t) = Y^p{t,n)h n {m). (12) 

n=0 

Function h n {m) equals to 

h n {m) = M [<5 m , a +6+...+d , (13) 

where 8 m n is the Kronecker delta. Due to the fact that S mn = — / e^ m ~ nS>!fi d(f, it 
is possible to write 



K{m) = —M 

Z7T 



e Hm-z 1 -b-...-tn)<P d(p = _ e imv Y[M[e-^ v ]dip. (14) 

7T . 27T J _ w 



Note that 



Therefore 



M [e-^] = ^ (e~ iLp + e iLp ) + 1 - a = a cos up + 1 - a. (15) 



h n (m) = —( e im ^(acos<p + l-a) n d<p. (16) 

Substituting this equation for h n (m), and earlier found expression for p(t,n) into 
the formula for f(m,t), we obtain 

°° +n -t/r i /-7T 

f(m, t) = J2 -^TTV- / e%mtf (« cos ^ + 1 - a) n rf^ = (17) 

£ — ' r n n\ L*K I - 
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27r ./-* 



n=0 

1 /-7T 

= J_ e -*/r f /" e ij7HP e t(acosvj+l-a)/T^ _|_ /" e -imp e t(a cos p+l-a)/-r^ 

= -e~' /r / cos(m^)e* (acoS¥,+1 " a)/T ^ = -e^" / cos(m^)e tocosv/r ^. 
n Jo tt J 

The last equation can be reproduced in a different form, using integral representation 
of Bessel function J m (z): 

J m (z) = / cos(m^)e i2COS ^. (18) 

Jo 

As a result, we obtain 

f( mi t) = i m e- at ' T J m ( K -i^j. (19) 
Taking into consideration that 

J m (iz) =i m I m {z) : (20) 

where 

00 r z \ 2k+m 

/ m ( 2 ) = yJ| — (21) 

^ A;! (Jfe + m)! 

are the modified Bessel functions, we write the final result for the probability of 
location of the particle at node m at time moment t: 



f(m,t) = (-l) m e-^I m f-^j = e ~^I m (f) . 



(22) 



Asymptotics of f(m,t) at t — > 00 with the asymptotic behavior of modified 
Bessel functions 

e z 

I m (z) = —== (1 + O {z~ 1 )) at ^^oo (23) 
y2nz 

considered, is of the form: 

/M = ^vS? {l + < r '» = Tssp f 1 + f 1 )) ' (24) 
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6. Thermal equations of the Universe evolution 

It is assumed that at high temperature (at early stage of the evolution of the 
Universe) it is possible to describe matter using field theory. Equation of state with 
good approximation corresponds to ideal quantum gas of massless particles (in the 
general case, it can be a mixture of ideal Bose- and Fermi-gases). In this theory, 
under the condition that the temperature T is far from mass threshold yet (radiatioro 
dominance, p = 3p), thermodynamic functions are given by the formulae: 



P = 3P = ^N(T)T\ (1) 



%N(T)T\ (2) 



where N(T) is the function related to the number of bosonic and fermionic degrees 

7 

of freedom (N(T) = A^(T) H — Nf(T)), p and p are equilibrium energy density and 

8 

pressure of the matter, s is specific entropy. All expressions are written in the unified 
system of units c = H — 1. 

Formulate the dynamic equations of the evolution of the Universe in terms of 
temperature. 

Note . The required equations are not Einstein's equations in the standard form 
(see Einstein equations, for example, in [31 H]). It is proposed to write the equations 
of the evolution of the Universe in Friedmann model using thermodynamic func- 
tions and temperature as function of time. It is possible to do this with the use of 
energy conservation law and condition of adiabatic expansion of the Universe in the 
framework of the standard cosmologic model. 

Instruction . Introduce auxiliary function 



<T) = ^, (3) 



where k — 0, ±1 respectively for flat, open and closed models of the Universe with 
time- dependent scale factor of a = a(t). 



2 By radiation here, we mean any relativistic object, including relativistic matter as well as 
photons. 
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Guth, Alan Harvey - American physicist 
and cosmologist who has first proposed 
the idea of cosmological inflation. In 2004, 
Guth together with Andrew Linde were 
awarded cosmological prize named after 
Peter Gruber for his work on the theory 
of inflation Universe. 



SOLUTION 

This problem of finding the time dependence of the temperature of the Universe 
was first formulated by A. Guth in [5|. To solve this problem, we would need to 
rewrite Einstein-Friedmann equation 



d\ 2 k 8n 



- +~ = ^Gp (4) 
a J a z 6 

in terms of temperature |2] (here G is the gravitation constant). It is considered 
that stress energy tensor of the Universe takes the form of energy momentum tensor 
for ideal liquid [3]. We shall use the energy conservation law 

> 3 H-4(« 3 ) (5) 

(the change in energy in a comoving volume element, d {pa 3 ), is equal to minus the 
pressure times the change in volume, — pd(a 3 )) and the fact that in the standard 
cosmological model it is supposed that the Universe undergoes adiabatic expansion 

= (6) 

(the entropy per comoving volume element remains constant). 

Let us write Einstein-Friedmann equation (j3J) in terms of temperature supposing 
that temperature value is far from the mass threshold (see e.g. |3]). We deal with 
matter which is found in thermodynamical equilibrium at almost all the time stages 
during cosmological expansion, so the chemical potential is considered to be zero. 

Taking into consideration ([T]) and auxiliary function (j3J let us represent eq. (j3J) 
in the form 

a \ 2 4-7T 3 

- +e(T)T 2 = —GN(T)T A . (7) 
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Taking into account equation of state p = 3p, from the energy conservation law 
we obtain the relationship 

a _ lp 

a""!? (8) 



(9) 



which in agreement with ([Tj) would take the form 

a f 1 N(T) 
a~~f~ AN(T) ' 

Using the condition of adiabatic expansion of the Universe we can find 

± = -~. (10) 

a 3 s 

Substituting into ffTOl) the expression for specific entropy (j2J), we obtain 

a f 1 N(T) 



a T 3N(T)' 



(11) 



Comparing (jH]) and ( II ip . we come to the conclusion that N(T) = 0. So the 
relation between the scaling factor and temperature should have the form 

- = ~t (12) 

Substituting this equation into eq. (jZJ), we obtain one of the dynamic equations of 
the evolution of the Universe: 

f\ 2 Air 3 

-J +e(T)T* = —GN(T)T\ (13) 

To write the second required equation, let us multiply both parts of the eq. fl2J) by a 3 , 

express — — — — ■ and substituting it into the auxiliary eq. ([3]), we obtain the equation 
a 2 T 2 



e(T) 



2tt 2 N(T) 
45 S 



2/3 



(14) 



where S = sa 3 is the total entropy in the volume defined by the radius of curvature 
a. We need to note that N, S hence e are constant in the considered temperature 
(or time) range (however a is not constant) due to the particle counting by allowed 
particles degrees of freedom thresholds. 

As the result, the dynamic equations of the evolution of the Universe in terms 
of temperature and entropy are the equations (TT3|) and (|14p . 
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7. Trapped electron 

Consider an isolated conducting sphere of radius R carrying the total charge Q. 
At the distance a > R from its center, there is a point charge q (qQ > 0). Find 
potential of the system tp(r) and the force F(a) acting on the point charge. Analyze 
the limit lim F(a), explain the obtained result. 

a->R+0 

SOLUTION 

The present problem can be solved using method of image charges. 

It is known that for any two point electric charges of opposite sign, it is always 
possible to find such a spherical surface that the resulting potential on it would 
be zero. Radius of the sphere and the distance from its center to the charges is 
determined uniquely if the values of the charges and the distance between them 
are known. So, the system under discussion (the point charge and the conducting 
sphere) is equivalent to the set of point charges. Thus let us place a charge q\ on a 
line connecting the center of the sphere O with the q charge at the distance of d in 
the direction of the charge q. One more charge, q = Q — g 1; we shall place at the O 
point, see fig. 1. 

Let us place the origin of the reference frame also in O point and direct OX axis 
to the point charge q (leftwards). Let us write the condition that the total potential 
(due to the charges q and gi) is zero at points (R, 0, 0) and (— R, 0, 0). 

Hence we obtain that d = R 2 /a, q\ = —qR/a, thus q = Q + qR/a. 

The total potential inside the full sphere equals qo/R and outside is given by the 
expression: 



q 

r 



qi qo 
n r " 



(1) 



The same in the Cartesian coordinates: 
q R 



y/(x - a) 2 + y 2 + z 2 a _ ^ 



+ 



- ^-) 2 + y 2 + z 2 



a/x 2 + y 2 + z 2 



(2) 



F(a)q / 






X 


< 


q, | 

V *T 
-\ > 







Fig. 1 
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The force on q acts along axis OX. Its projection is given by 



/ 



qQ , q 2 R 



a 2 ; / 



(3) 



It is easy to check that if we represent a as R + A, in the limit at A — > 
the expression for the force would become F ~ — q 2 / (2A) 2 which corresponds to the 
interaction force of a point charge with non-charged conducting plane. 
qQ 

At a — > oo, -F — > — — (Coulomb force, as expected). 
a 2 

At a — > R, no matter how large Q is, even if Qq > 0, we have F — > — oo, 
i. e. the force becomes very strong and attractive! 

Let us find the distance at which the effect of attraction starts. Let us write the 
interaction force ([3]) in a dimensionless form: 



a 



m = -5 - 



2s 2 - 1 
sHs 2 - 1) 



(4) 



where 



« = _>0, 8 = — >l, f(s) = —F x (a). 
q R q 2 



(5) 



The force occurs to be zero at a distance a = SqR, where s > 1 satisfying the 
equation 



2s 2 - 1 
s(s 2 - 1) 



a. 



(6) 



Function f(s) graphs are shown in fig. 2 at a = {2; 1;0.5}. In these cases sq = 
= {1.43; 1.62; 1.88}, and force f(s) reaches its maximum {0.43,0.15,0.05} when 
s = s max = {1.79,2.07,2.46}. 

This explains why electrons can't escape from metals, even though they are 
repelled by the other electrons. If an electron manages to escape and gets to a small 
distance from the surface of metal, the other electrons conspire to bring it back by 
rearranging themselves in such a way to create a huge image charge which attracts 
the electron back to the metal with a strong electric force! 



Mathematical Physics: Problems and Solutions 



25 




Fig. 2 



8. X-sector @ 

Consider the following model of Higgs sector with two doublet scalar fields <pi 
and 02 transformable as SU(2) doublets, with the weak hypercharge generator of 
Yy/ = 1, and with each component of the doublet being a complex scalar field. 
Suppose both fields acquire parallel vacuum averages (vevs) of the type 

<*> = 7iU) (i = 1 ' 2) ' (1) 

with the values Vi, v 2 (these vacuum averages lead to gauge bosons mass matrix as 
in the Standard Model with the replacement Q v 2 = vf + v%). The most general form 
of the potential energy function (potential) for a model with two Higgs doublets is 
rather complicated!!. However, the model hermitian potential possessing the main 
properties can be written in the following form: 

v{<f>xM = -»M<h.) - vMfa) - fiz^M - (a4)* (^0i) + 

+A 1 (<^ 2 + A 2 (0 2 2 ) 2 + A 3 (^ 

where [jl\ 2 an( i ^5 ma y be the complex numbers. 
3 Some extension of Peskin & Schroeder [B] problem 20.5. 

4 Recall that in the Standard Model (SM), the gauge bosons masses come from the term \Dcf>\ 2 
in the Lagrangian, where we set <p equal to its vacuum expectation value v. 

5 It is only for reasons of simplicity that the SM contains just a single Higgs doublet. Supersym- 
metric extensions of the SM typically contain two or more Higgs doublets, and singlets. 
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(a) Obtain the conditions that for the direction in field space at given configura- 
tion of vevs (Q]) the potential is bounded below at large field values. (The analogue 
of A > for the theory with a single Higgs doublet.) 

(b) Find the conditions to impose on the parameters \i and A, so that the config- 
uration of vacuum averages (JT|) gives strictly local (locally stable) minimum of this 
potential. 

(c) In the unitary gauge (rotation to canonical form), one linear combination of 
the upper components 0i and 02 nulls, while another becomes a physical field. Show 
that charged physical Higgs field is of the form: 

H + = 0+sin/3-0+cos/3, (2) 

where is defined by the relation 

tg/3 = -- (3) 

Vl 

(d) Investigate if the CP invariance breaks in the given potential. Substantiate 
the obtained results. 



SOLUTION 

(a) First of all, note that fif =1 2 and Aj=i...4 are all real due to the fact that 
Lagrangian is Hermitian. But [i\ 2 and A5 may be the complex numbers. 

Second, find the condition under which both <fti and 2 have parallel non-zero 
vevs (recall that in the case of a single Higgs doublet the gauge symmetry allows 
the vacuum expectation value (vev) to be taken in the form (0, v), with v real). Now 

Wo 



use a rotation to make (0i) = —= ( , where v\ is real. Note that after we have 

V2\vi J' 

done this we have used all of our gauge rotation freedom, so the vevs of 2 are still 
completely general, i. e. (02 ) = — ^ ( V Ji V So we rewrite the potential in terms of 
the two (complex) vevs of 02, and of the V\. Letting v\ = v' 2 2 + v 2 2 , we find 

t t jp/ \ Ml2 11 (^12) //* 1 ^4 2/ //* ll\ 1 "^5 2 //2 . ^5 2 li* 2 ( a\ 
V = F(V 1 ,V 2 ) - — V X V 2 — V ^ V 2 + -J V l( V 2 V 2 ) + -^ V 1 V 2 +-g-<^2 > ( 4 ) 

where we denote F(vi,v 2 ) the function dependent on values V\ and v 2 only. 

First we want to answer two questions: (i) What is the condition that v 2 is real? 
(ii) What is the condition that v' 2 is zero? If we enforce these conditions, it ensures 
that the vevs take the form of (fTj) , where both v\ and v 2 are real. This is what it 
means for the vevs to be 'aligned'. So, how to find these conditions? It becomes clear 
if we rewrite v 2 = ae td . The potential written in terms of a and 9 is given by: 

V(a, 9) = (Stuff not dependent on 9 or a) — Re/i 12 t>ia cos 9 + lmfi 12 Viasm9 + 

+ -\<±v 2 a 2 + -ReA 5 t; 2 a 2 cos 29 — -lm\ 5 v \a 2 sin 29. (5) 
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The reality of v 2 is ensured by forcing 9 = 0, it. We can see that 9 = or it will 
be a stable minimum of the potential if the second derivative of expression §5§ is 
positive at these values of phase 9, i. e. 

ReA 5 - Reu? 2 — < 0, (6) 
via 

and in the limit of large field values 

ReA 5 < 0. (7) 

In this case the Eqn. ([5]) will be minimized for cos 29 = 1. We must show that v' 2 = 0. 
Recall that the sum v' 2 + v 2 2 is fixed to be v\. So if we can arrange the potential 
so that it is energetically advantageous for all the vevs to go into v 2 2 , we are done. 
This is equivalent to saying that we want V(a, 9) to be minimized for a — > oo. This 
is accomplished if 

A 4 + A 5 < 0. (8) 

So, together, eqns. ([7j), ([8]) guarantee that the vacuum expectation values can be 
aligned. 

(b) Now what is still required, is to show that this is a stable minimum. Using 
aligned forms for the vevs, rewrite the potential in terms of v\ and v 2 . What addi- 
tional conditions on the parameters are necessary to guarantee a stable minimum? 
Stability is equivalent to saying that there is a positive mass squared for fluctuations 
about the minimum. In other words, we examine the mass matrix Hessian 



H 



( d 2 V(v 1 ,v 2 ) d 2 V{v 1 ,v 2 ) \ 

dv\ dv 1 dv 2 
d 2 V(v u v 2 ) d 2 V( Vl ,v 2 ) 



\ dv 2 dvi dv 2 ) 



in minimum 



In order to have a stable minimum the matrix of second derivatives needs to be 

positively definite. This does not mean that all second derivatives need to be positive. 

Both eigenvalues, e^, need to be positive to ensure the stability of the minimum. 

The above matrix is to be evaluated at the minimum, setting to zero the derivatives 

dV dV 

evaluated at the vevs (i.e. where — — = — — = 0). Since invariant Ty(H) = e\ + e 2 

OVi ov 2 

and Det(if) = e\e 2 , we can simply require 



d 2 V( Vl ,v 2 ) | d 2 V{v u v 2 ) 
dv\ dv 2 

Straight-forward algebra gives: 



> 0, and Bet(H) ^ 0. 



2A!^ + Reul 2 — - Re/i^ 2 + A345W1W2 

H = I Vl 2 Vl ] ■ (10) 

- Re/i? 9 + \345V1V2 2\ 2 v\ + Re/i 12 — 

v 2 
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We denote A 345 = A 3 + A 4 + ReA 5 . 

This shows that stability conditions are equivalent to: 



2X 1 v'i + 2X 2 vl+ ( — + — ) Re/4 > 0, (11) 



and 



4AiA 2 + 2Ai^Re/4j + 2A 2 ^|Re/^ 2 ^ (A 345 ) 2 - 2A 345 Re/i 2 2 , (12) 

by taking first derivatives with respect to all the scalar fields, and setting them equal 
to zero (i.e. in minimum). 

In particular case Re/i^ = we see simply 

4AiA 2 > (A 345 ) 2 , and Ai > 0, A 2 > 0. (13) 

In addition, we will requite concavity, which is implied to be positive, which 
assures us that we're at a minimum and not at a saddle point. The tricky part of 
the problem is to show that indeed it is possible to have the vevs parallel, i. e. of 
the form of (CE]).E] To do this, the best way is to use the SU(2) rotation to force the 
vev of 0i to have the right form. Then all we need to show is that the potential 
is minimized (i. e. the appropriate derivatives satisfy the conditions stated above) 
when <p2 takes the right form. 

Thus the conditions for a local stable minimum of the potential in this problem 
are ©, ©, (HU), (ED, see also 0. 

(c) In the SM, there are 4 degrees of freedom in the Higgs doublet, three of 
which are consumed by the W + , W~ and Z°. When there are 2 Higgs doublets, 
they contain in total 8 degrees of freedom, so the 5 remaining after goldstones are 
consumed. Now, the vev given breaks the SU(2) x U(1)y symmetry to U(l) em , 
therefore in general there would be 3 goldstone bosons, and 5 physical Higgs fields. 
Two of these remaining degrees of freedom are charged, and three are neutral. The 
task in this problem is to determine which two charged degrees of freedom are eaten, 
and which two charged degrees of freedom remain. The simplest way to do this is 
to consider the components of the two Higgs doublets as being part of a larger 
vector. An orthogonal transformation will rotate the different components amongst 
themselves. In particular, we can find the basis where the vev is entirely in one 
neutral component. The charged piece associated with this neutral component is 
the would-be goldstone boson that is eaten. The charged Higgs is the piece which 
is orthogonal to this goldstone boson that is eaten. So, considering the neutral 
components (which have the vevs), we have: 

V i \ _ ( cos @ ^f 3 
i> 2 J I — sin f3 cos J I 



(14) 



6 Thus, in fact, condition ((6]) is equivalent to positive sign of squared mass for third Higgs 
(pseudoscalar) boson m\ = — ufReAs+Re/^ — > an d condition at large field values is equivalent 

v 2 v 
to positive sign of squared mass for charge Higgs boson m 2 H ± = — ^"(A4 + R- e -^5) + Re//f 2 — ■ 
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So, putting all of the vacuum in one Higgs field0 (say 0^), we get 

/ cos/3 sin/3 \ / v x \ ( v \ . , 

V - sin/3 cos/3 J \v 2 ) \ )' 1 ' 

It follows from (fT51) that tg/3 = — . So, the whole vev lives in the <p'i field, which 

Vi 

means that the charged component of 0'i is the Goldstone Boson eaten by the W. 
That means <f)[ = 0i cos /3 + 02 sin /3 is the Goldstone, and 

0' 2 = H + = 0+sin/3 -0+cos/3 (16) 

is the physical charged Higgs field. 

(d) In the models with two doublets of scalar fields CP invariance can be vio- 
lated by the terms of the potential containing (0^02) or (020i) with the complex 
parameters //f 2 and A5. In the case of real parameters CP invariance is not broken. 



9. By the cradle of LHC | 

Generally, it is possible to describe a scattering experiment in the following way 
(see fig. 3): 

1) sufficiently wide uniform beam of particles is prepared so 
that it is possible to assume the momentum of each particle be 
equal to po = hk, where k is wave vector, h - Planck constant; 

2) this beam of particles is directed to stationary target con- 
sisting of identical particles; 

3) at certain distance from the target, products of reaction 
of the particles from the beam with the particles forming the Fig. 3 
target, are registered at different angles. 

Thus at sufficiently large distances from the target the wave function of the 
particles is the superposition of plane incident wave -0m = exp{ik r} and spherical 
scattered wave ip = exp{ikr}/r. Here, r = xi + y] + zk, r = \ f\. 

h 

(a) Calculate flux density of the probability j = (^*V^ — ipVip*) for the 

2mi 

wave functions ^ = e lh ' x and ip . Here, m is the mass of the particle, and x = ix. 

(b) Illustrate the obtained result with the help of the graph: draw the pattern of 
the vector field j (lines of the j vector) in both cases. To build such graph, use any 
available computer software suitable for building graphics (plots). 

(c) Prove that V • = 0, V • j ~ 5{f). 



7 Once one goes to the proper basis for describing the Higgs mechanism, there is really only one 
doublet (in the above case it is (/)[) that acts as the Higgs and has three of four degrees of freedom 
that are eaten. So, we can say that it is a model with "two-complex scalars". 

8 Large Hadron Collider. 
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SOLUTION 

(a) 

JW = ^(VCWw - v-v^: 



2mi 



i e 



m 



h 



J 



h 



W> vv* ) 

O 2mW © O O © 



n. 



r dr r 



dr r 



hk n, 
m r' 



2mi 
Here, n r = f/r. 

(b) 

The field in both cases is represented by a graph in x—y plane at z = (in 
Mathematica realization, see fig. 4). 



In[1]:= Needs [ "VectorFieldPlots ] ; 

ln[2]:. VeotorFieldPlot[{l, 0}, {x, -1, 1}, {y, -1, 1}] 



ln[3]:= VectorFieldPlot [ - 



{x, y> 



{X, -1, 1}, {y, -1, 1}, 



(* 2 + y 2 ) 

ScaleFunction -» 

WWW \ \ t f / / / // 

ww \\ \ t f / / / / / 

WW \ \ \ / / / / / /s 

v.v.n\\ \ \ / / / / sss 

^^^w\\ \ / / //^^^ 

^*^-*W\ \ I / / y^^^^ 



s/ / i \ \w^-^— 

/ / i \ \ \ w^^ 

"// / / / \ \ \\xw 

'// / / / / \ \ \ \\xw 

// / / / / / M \ \ 

/ / / j I I I \ www 



Fig. 4 

It is possible to see that the spherical wave ib is diverging from the origin of 

O 

coordinates. 



Vector field is uniform, = i 



const, so that V • = 0. 



It is possible to represent vector j in the form: 



2 

O m r z 



hk n r hk f hk^l 



m r- 



m r 



Hence 



o 



hk A 1 hk , . , . 
— A- = — 4vr5(f). 
m r m 
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10. 'Whipping Top-Toy' from Samara § 



Rigid ball of the mass m with radius R rests on smooth rigid horizontal surface. 
Center of mass C of the ball is at the distance of / from its geometric center O. Mass of 
the material is symmetrically distributed along the volume of the ball relatively to 
OC axis, and also any plane containing that axis. Moments of inertia of the ball 
relatively to OC axis and any axis passing the center of mass and perpendicular to 
OC axis, are equal respectively to Jo and J. During certain time period, the ball 
is accelerated around the static vertical axis passing its center O. The moment of 
time when the action of the "accelerating" forces is finished, is chosen as the time 
origin. At this moment the ball has the angular velocity uto, directed vertically up, 
and OC axis makes some angle e with the rotation axis of the ball. The center of 
mass of the ball is lower than its mechanical center (see fig. 5) . 




Fig. 5. Rotating ball at time moment t = 



Study the motion of the ball at t > 0. Obtain the equation of motion of the ball 
and find the integrals of motion (conserved values). Show that the center of mass 
of the ball is lifting up (at a certain relationship between the parameters), and find 
out, whether it can approach its top position at which the OC vector is directed 
vertically up. 

SOLUTION 

In order to describe the motion of the ball we introduce an inertial frame of 
reference 5" (with axes X, Y, Z), that is at rest relatively to the horizontal surface, 
and an noninertial frame S' (with axes X', Y', Z'), that is rigidly bound to the ball. 
The origin of the frame S coincides with the initial position of the center of the ball, 
the axis Z is directed vertically and the plane YOZ is chosen so that it contains the 
axis OC at the initial moment t = 0. The origin of the frame S' coincides with the 
mass center of the ball and the axis Z 1 - with the axis OC, so that at the moment 
t — the axis Z' makes an angle e with the axis Z. Additionally, the axis X' at this 



9 Place in Russia where this competition is held and assessed. 
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moment has the direction similar to that of the horizontal axis X, and the axis Y' 
belongs to the plane YOZ (see fig. 6). 



z 



ft- J 

z'\ V 

\ E 


"\ t= 




\— — >y 






7 


i 





Fig. 6. Coordinate axes and the origins of S, S' at the moment t = 

The reference frames S and 5" being chosen, the Euler angles tp, 9, ip that define 
the orientation of the ball (and the frame S') relatively to the frame S, are given by 
the following values at the initial moment: 

y>(0) = 0, 0(0) =e, ^(0) = 0. (1) 

Components of the angular velocity Q of the ball in the frames S and S' at any 
instant of time are given by the cinematic Euler's formulas: 

u x = 9 cos tp + ip sin 9 sin tp, 

u y = 9 sin tp — ip sin 9 cosy, (2) 
u z = tp + ip cos 9, 



bj x i = tp sin 9 sin ip + 9 cos "0, 

Uy' = ip sin 9 cos ip — 9 sin ip, (3) 
u z > = tp cos 9 + ip. 

For the initial values of the components of the angular velocity in S the equa- 
tions §Z§, with the account of (CQ), give 

w*(0) = 0(0), 

Wtf (0) = -^(0)sine, (4) 
u t (0) = (p(0)+ip(0) cose. 

On the other hand, it is known that the initial angular velocity uj(0) = uj has 
the same direction as the axis Z and, hence, 



uj x (0) = u y (0) = 0, u z (0) = uj . 



(5) 
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Comparing (j3J) and (jSJ), one finds the initial values for the time derivatives of the 
Euler angles: 

0(0)= cu , 0(0) =0, <H0) = 0. (6) 

As seen from the pictures, the cartesian coordinates x m , y m , z m of the mass 
center of the ball in S at the initial moment of time have the following values: 

x m (0) = 0, y m (0) = / sine, z m (0) = -I cose. (7) 

For t < (until the rotating forces action has stopped) the center of mass of the 
ball moves in the horizontal plane along the circle of the radius / sin e with the axis 
Z passing through its center. Therefore, the velocity v rn of the mass center of the 
ball in the frame S at the initial instant of time is given by 

v m (0) = p f m (0)], (8) 

where f m (0) is the initial value of the radius- vector r of the mass center of the ball 
in S (that is the vector with the components x m (0), y m (0), z m (0)). From (j5J), (J7|) 
and (jHJ) it follows that the initial velocity of the mass center of the ball in S is 
directed opposite to the axis X and is equal to 

v m (0) = u lsine, (9) 

or 

im(0) = -u l sine, y rn (0) = z m (0) = 0. (10) 

The equations of motion of the ball for t > can be derived as the Lagrange 
equations of the 2-nd kind. The ball is a mechanical system with five degrees of 
freedom. One can use the Euler angles tp, 8, ip and coordinates x m , y m of the mass 
center in S as the independent generalized coordinates of the ball, the coordinate 
z rn at any instant of time being defined by 

z m = —I cos 9. (11) 

The Lagrange function has the form 

L = T-U, (12) 

where T is the kinetic energy of the ball and U is its gravitational potential energy 
(these values are defined in the frame S). 

The kinetic energy of the ball may be expressed in the form 



(13) 
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where r is its kinetic energy of rotation. As the axes x', y', z' are directed along the 
principal axes of inertia of the ball, then it follows 

r = -{J x >u 2 x , + J y ,uj 2 y , + J z >u 2 z ,), (14) 

where J x >, J y >, J z i are the moments of inertia of the ball relatively to the axes x', y', 
z', with, according to the problem statement, 

Jx> = Jy' = J, Jz> = Jo- (15) 

Inserting the expressions Q and (TT5I) into (TT4)) one gets 

t='U^ 2 sin 2 9 + 9 2 ) + ^{<p cos 9 + ^) 2 . (16) 

The expression for the velocity's square of the mass center of the ball with the 
account of (ITT]) gets the form 

v 2 m = x 2 m + y 2 m + l 2 sm 2 99 2 . (17) 

From (H3J, (HSJ and (HZJ it follows that 

T = j (x 2 m + ijl) + jl 2 sin 2 99 2 + J -{v 2 sin 2 9 + 9 2 ) + 

+ y ^ 2 cos 2 ^ + 2^cos^ + ^ 2 ) . (18) 

The potential energy of the ball with the account of ( TTTj) is given by the expression 

U = -mgl cos 9. (19) 

Inserting the equations (|18|) and f[T9"j) into (112)) one gets the following explicit 
expression for the Lagrange function of the ball: 

777 

L = —(x 2 m + y 2 J + -(ml 2 sin 2 9 + J)9 2 + -( J sin 2 + J cos 2 6# 2 + 

+ + Jo cos 9ipip + mgl cos 9. (20) 

It is known that the Lagrange equations of the 2-nd kind for a system with ideal 
holonomic constraints and S degrees of freedom under the absence of dissipative 
forces have the form 

KS-lr ^" (21) 
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where qj(j = 1, 2, S) are independent generalized coordinates of the system. 
If qj is a cyclic coordinate of the system (that is, the Lagrange function doesn't 

dL 

depend on it), then the corresponding generalized momentum = — — is the inte- 

d qj 

gral of motion. 

According to (12 Op . x m , y m , tp, ip are cyclic coordinates of the ball. Therefore, the 
generalized momentums 



dL 

dx m 

dL 



(22) 



my m , (23) 



dL 

dip 



J sin 2 6^0 + J o cos6»(cos6»0 + tp), (24) 



dL 

-j = J (cos9<p + i>) (25) 
dip 

are the integrals of motion of the ball that are constant and are defined by the initial 
conditions ([I]), (jHl), © an d ffTOj) . Writing down the corresponding conservation laws, 
one gets the following system of the 1-st order differential equations: 

x m = — uqI sine, (26) 
Vm = 0, (27) 
J ip sin 2 6 + J u cos e cos 6 = u (J sin 2 e + J cos 2 e) , (28) 

(p cos 6 + ip = ujq cos e (29) 

(when deriving the equation (|28j) . that presents the conservation law for the gener- 

dL 

alized momentum = ——, the equation (1290 was used). 

dip 

From (12"6"|) and (12T)1 and the initial conditions (J7j) it follows that 

x m = — UqU sine, (30) 

y m = /sine. (31) 
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Hence, the mass center of the ball moves in the vertical plane that is normal 
to the axis y and that passes through the initial position of the mass center. The 
horizontal component of the ball's velocity stays constant and equal to the initial 
velocity v m (0). 

Note that the equations (126]) - (128]) may be also derived from the conservation 
laws for the components of momentum and angular momentum of the mechanical 
system in an inertial frame. In order to do this one has to consider the directions 
of external forces, applied to the ball. These forces are the reaction force N that is 
directed vertically up and the gravitational forces that are applied to every element 
of the ball and are directed vertically down, their resultant force mg being applied 
to the mass center of the ball. 

As the sum of the external forces that are applied to the ball, is directed vertically, 
the X and Y components of the ball's momentum 

P = mv m (32) 

are conserved for t > 0. The corresponding conservation laws, with the account of 
the initial conditions (|10j) . give the equations (12"6"|) and (|27j) . 

The angular momentum M of the ball in the frame S may be presented in the 
form 

M = [fP] + jx, (33) 

where fl is the intrinsic angular momentum momentum of the ball (that is the angu- 
lar momentum of the ball in the mass center reference frame that moves translatory 
to S). The components of ft in 5" are defined by 

/V = Ju x >, [My/ = Joy, /v = JqUs z >. (34) 

As the moments of all external forces in S that are acting on the ball are directed 
horizontally for t > 0, the projection M z of the angular momentum on the axis 
Z is the integral of motion. Taking into account ( 130]) . (I3ip . ([TT]) . (126]) and ( 127]) . it 
easy to check that Z-component of the vector [fP] has the constant value equal to 
ma;oZ 2 sin 2 £. Hence, with the account of (133]) . it follows that 

M z = mcuol 2 sin 2 e + fx z , (35) 

or that Z-component of the ball's rotate momentum \x z is an integral of motion. 

Using the well-known linear transformation law for the components of an arbi- 
trary vector for a frame rotation and expressing all the coefficients of such transfor- 
mation as the functions of Euler angles, one may show that 

ix z = sin 9 sin i/;[x x i + sin 6 cos ijjLiy> + cos 8fx z i . (36) 

Inserting the equations (154]) . where the angular velocity components are given 
by the formulas ([3]), into (I36p . one gets 



jx z = (J sin 2 9 + J cos 2 9)Lp + J cos 9ip. 



(37) 
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Comparing the equations (I37p and (I24p . one sees that \i z coincides with the gen- 
eralized momentum p v . Hence the conservation law for fi z , written with the account 
of the initial conditions ([1]) and (Q and of the equation f l29|) . is the equation fl28|) . 

Inserting the equation (120]) into (l2~Tj) and letting qj = 0, the Lagrange equation 
for the ball that corresponds to the generalized coordinate is obtained: 

(ml 2 sin 2 9 + J)9 + ml 2 sin cos 00 2 + ( J - J) sin cos 0yi 2 + (38) 

+Jq sin 0<^-i/? + mg/ sin 6 = 0. 

Note that the equations (}2"8"j) and (T2"9"j) form a system of linear algebraic equations 
relatively to (p and ip with the coefficients that depend only on the angle 0. Having 
solved this system, one gets the dependencies of these generalized velocities on 0: 

(f = ( J sin 2 e + J cos 2 e — J cose: cos 9), (39) 

J sin # 

ip = — ^^[cos^Jsin 2 ^ + J cos 2 9) - cos #( J sin 2 e + J cos 2 e)]. (40) 
J sin 

After having inserted the expressions (1391) and fj40|) into (138|) one may get the 
nonlinear 2-nd order differential equation that defines the dependency 9{t). In order 
to get this dependency one may also use the conservation law for the total mechanical 
energy of the ball 

E = T + U. (41) 

This quantity coincides with the generalized energy of the ball and is conserved due 
to the absence of both dissipative forces and explicit dependency of the Lagrange 
function on time. 

Setting the explicit expression for the energy of the ball that follows from f fl8|) . 
( TT9"]) and fl4"T|) . equal to its initial value, that is found from ([T]), §6§ and f lTU]) . and 
taking into account the conservation laws ( l2T)|) and ( l2~T|) . one gets 

(ml 2 sin 2 9 + J)9 2 + (J sin 2 9 + J cos 2 9)<p 2 + J ip 2 + 2 J cos 9Lpip - 2mgl cos 9 = 

= (J sin 2 e + J cos 2 e)uq — 2mgl cos e. (42) 

Substituting the expressions (13"9"|) and fl4"0l into the last equation, one obtains 
the following nonlinear 1-st order differential equation that defines the dependency 
9{t): 

sin 2 fl(l + /3sin 2 fl)^ 2 = oo 2 (cose -cos 0)0(0), (43) 
where 0(0) is quadratic in cos0 and has the form: 



0(0) = a + 2oi cos 9 + a 2 cos 2 0, 



(44) 
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where do — d 2 are given by: 

a = -2/?7 - - ((a - l) 2 cos3e + (a + l)(3a - 1) cose) , (45) 

ai = - ((a 2 - l)cos2e + a 2 + l) , (46) 

a 2 = 207, (47) 

Dimensionless parameters a, /3, 7 in equations (145]) - (J47J) and equation (143]) . are 
related to the known parameters in the following way: 

o = y, (48) 
= ^, (49) 

(here u m = y/g/l is the cyclic frequency of a flat simple pendulum of length I). 

The region of motion of the ball top is defined by the sign and the roots of 
the function 0(6*). In order to investigate them one has to specify the range of the 
parameters that are present in the function. As J, Jq > and as for a particular case 
of the ball being a rotator the moments of inertia have the values Jo = 0, J 7^ 0, we 
conclude that a is not negative. On the other hand, it is known that none of the 
principle moments of inertia is bigger than the sum of the other two, which gives 
J < 2J. Therefore, we conclude that < a < 2. Further, as it is easy to notice, 
< < 00 and < 7 < 00. 

Taking into account the equations ( 145]) - (|47|) . one finds that: 

< 01 < 2, (51) 



< o 2 < 00. (52) 

Concerning do, one may notice that this parameter is not bound from below and 
approaches its maximum value with respect to /?7 at /?7 = 0. Further analysis shows 
that for < a < 2, < s < ir the parameter do reaches its maximum value, that is 
equal to do = 4, at a = 2, e = ir. For small e the parameter Oo is maximized with 
respect to a at a = and increases monotonically with e. Hence, we finally conclude 
that 



— 00 < a < 4. 



(53) 
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From the form of the equation (j43p one can conclude that for 0(e) > the 
motion of the top will always go on in the area where cose > cos 9 or, that is, where 
9 > e. Correspondingly, for 0(e) < the mass center will always be lower than in 
its initial position during the motion. Finally, for 0(e) = the motion will be going 
on at constant 9, that is, the mass center of the top will be moving at constant 
height (however, this regime is unstable with respect to small variations in a,/3,j). 
The last result follows from the fact that for 0(e) = 0, as we shall see further, the 
right-hand side of ( 143]) is negative everywhere except for at 9 = e. 

From the equations (j44|) and (T4"5"]) - (147|) it follows: 

0(e) = -2sin 2 e((a- l)cose + /?7) (54) 

Hence, the condition that the mass center lifts up during the initial stages of 
motion is: 

0(e) > (a - l)cose + /?7 < (55) 
Note, in particular, that for the tops with a > 1 the mass center will always 

7T 

be moving downwards under the condition e < — irrespectively to the kinematic 
parameters of the problem. 

Remark . Note that in [2] the case for small e is considered in details and analytical 
solutions are obtained, and the analysis of trajectories is carried out also. 

Let us now study the behavior of the roots of the parabola, defined by the 
dependency 0(cos#). As ci2 > 0, the parabola is open upwards (the case CI2 = will 
be understood as the limiting one). The minimum of the parabola corresponds to 

cos 9 equal to z m = — — . As from the equations (I45p - (I47p , (l48p - (}5T)j) it follows 

that di, 02 contain independent parameters, one concludes that 

— 00 < z m < 0. (56) 

Further, from (1441) one gets that: 

2 

m = 0(z m ) = + o . (57) 

02 

Substitution of (j45]) - fT4TP leads to the expression: 

0m = -2/3 7 - ( a2c ° s2 £ + Sin2£ ) 2 _ CQS5(a 2 _ gin 2 5(1 _ a) 2y (5g) 

8/37 

As the parameters a, /3j may be considered as independent, one can say that 

m < —(a 2 cos 2 e + sin 2 e) — cose(a 2 — sin 2 e(l — a) 2 ). (59) 

The last expression becomes the equality when (3~f = -(a 2 cos 2 e + sin 2 e). After 
elementary transformations one arrives at 

6 m < -2cos 2 ^(l + (a-l)cose) 2 . (60) 
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Therefore, it is seen that the minimum of the parabola defined by the dependency 
0(cos#) is a strictly nonpositive quantity that achieves its maximum value equal 

to zero, under e = ir (and, simultaneously, f3j = ^« 2 )- Hence, the binomial 0(cos6 l ) 

does always have two real roots with respect to cos 9. Let us introduce the notation 
z + , z_ for the bigger and the smaller root correspondingly: 



Further, notice that: 
d6(cosi 



dcos^ 



* = ± - * (61) 

02 V fl 2 02 



4/37cose + a 2 cos 2 e + sin 2 e. (62) 



7T 

The last equation is strictly nonnegative under e < — . Hence, for a motion that 

starts from the positional where e < — , one can state that the region of motion will 
be defined by 

max[— l,z+] < cos 9 < cose (63) 
when the condition (1551) is satisfied and 

cose < cos 9 < min[^ + , 1] (64) 

when it is not satisfied. 

Considering the equations fl63|) - (|64|) and the equation (I61|l . it is easy to get the 
condition that at its highest position the mass center belongs to the same horizontal 
plane as the center of the ball. In order for this to happen the condition (155 p and 
the equation z + = must hold, which is equivalent to do = or 

8/?7+ ((a- l) 2 cos3e + (a + l)(3a- l)cose) =0. (65) 

Note that the last condition can be satisfied only if it admits real roots for 
a (the other parameters being fixed) and only if one of the roots is in the range 
< a < 1. The roots for a may be transformed to the form: 



ai 2 = — tt~ (— sin e ± \ sin e — 2/3^ cose) (66) 
cos 2 e v 

The condition of the existence of two real roots implies the inequality 

sin 2 e-2/3 7 cose > (67) 



7T 

10 Further in the text, the condition e < — will always be considered true. 
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The condition a > implies than only the bigger root may be a physically 
meaningful value and that the following inequality holds: 

sin 2 e cos 2 e — 2/3^ cose > 0, (68) 



Note that it already includes the inequality (1671) . The condition a < 1 leads to the 
inequality 

cos 2 e + 2/3 7 cose > 0, (69) 

which is satisfied automatically. 

We can finally conclude that for a top that was initially at the position where e < 

7T 7T 

— , the highest position for its mass center will be 9 = — if the conditions (16"6"|) . (I55p 

hold, the necessary condition for the first one being (16"5]) . 

Let us now study the values of the root z + . When (j55p holds (the top lifts up), 
it is always true that z + < cose. The lower bound for z + may be derived after 
inserting the equations (l4"5j) - fT47j) into (1611) . which gives, after transformations: 



where 

A 2 = sin 2 e + a 2 cos 2 e, (71) 
r = cose(a 2 cos 2 e + (—1 + 2a) sin 2 e), (72) 



fi = 4/3 7 . (73) 

As in (170|) the parameter is independent with respect to the other ones, one 
can minimize this expression with the account of fi > 0. It is easy to show that the 
expression for z + is minimal in the limit \i — > 0, converging to the value 

T 

2 ; +min C^O 

Now, considering a as a variable, it is easy to see that the equation z +min is 
minimal for a = and in this case is equivalent to z +min \ a=Q = — cos e. Therefore, 
we conclude that, irrespectively to the parameters of the problem, if the mass center 
lifts up at the beginning of motion, it is always in the range [e, ir — e]. 

Similarly, for the downwards motion the minimal value for z + is cose. The maxi- 
mum value may be derived from the equation (1701) . Maximizing this expression with 
respect to /i, one obtains that z + reaches its maximum value that is equal to unity 
in the limit fi — > oo. Therefore, one can conclude that for any initial position of the 
top the parameters of the problem may be chosen in such a way that the top may ap- 
proach arbitrarily close to the position 9 = 0. This result, in particular, corresponds 
to a simple case u$ = 0. 

In conclusion, notice that the equation (l43p may be easily integrated in elemen- 
tary functions by the use of the variable change r = cos 9. 
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11. Laplacian A spectrum on a doughnut 




Fig. 7 



Consider a torus made of a rectangular block 

< x < a, <y <b 

with glued opposite sides (see fig. 7) where the 
identical arrows mark the sides to be glued to- 
gether. For sufficiently big ratio a/b, it is pos- 
sible to implement such torus nearly without 
deformations, as the surface of a doughnut in 
3D space. 



If inside the doughnut magnetic field with the flux $ x is cre- 
ated, and also magnetic flux $2 is passing through the hole of 
the doughnut (see fig. 8), then the wave function of the station- 
ary state of the charged particle with the charge e and mass 
m, on the surface of the torus, is the eigenfunction y) of 
the operator 




Fig. 8 



2m 







dx ch 



i A x (x,y)) +[—-i—Ay(x,y) 







dy ch 



(1) 



where A(x,y) is a vector potential, h is Planck constant, and c is speed of light. 
Function ty(x,y) satisfies to periodic boundary conditions: 



*(x,0) = *( Xl b), *(0,y) = y(a,y), 



(2) 



*'Jx,0) = *'Jx,b), %(p,y) = *' x (a,y). 



(3) 



Suppose the magnetic field turns to zero on the surface, we can make cuts of 
the surface and, using gauge transformation, nullify A on the surface of the torus. 
Then the wave function would be discontinuous at the cuts (such cuts are possible 
to make exactly on the border of the rectangle). It is also clear that \^(x, y)\ 2 should 
not change. 

After the gauge transformation, the task to find stationary states is modified and 
requires to find eigenfunctions ip(x,y) of the operator 



H x = - 



1m 



d_ 

dx 



d_ 

dy 



(4) 



with phase shift (lagging) boundary conditions: 



ip(x,b), ip(0,y) 



(5) 
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(6) 



(a) Find eigenfunctions and eigenvalues of the operator H\. 

(b) What is the relationship between the fluxes $1, $2 and phase displacements 



SOLUTION 

Let's make mentioned above cuts of doughnut surface and nullify vector poten- 
tial A' with the help of gauge transformation 



A'(x,y) = A-Vf(x,y)=0. 



(7) 



Thus the Hamiltonian Hq is simplified. In spite of the fact that the wave function 
would undergo certain changes, it is clear that the probability \^f(x, y)\ 2 of particle 
location should not be changed by gauge (gradient) transformation, i. e. the new 
wave function ip(x,y) differs from the old one ty(x,y) only in phase factor [8]: 



ip(x, y) = V(x, y) exp ( - z^/fo y) 



1 e 



(8) 



where f(x,y) is the function of gauge transformation, ip(x,y) is the eigenfunction 
of the new modified Hamiltonian Hi, that is in fact Laplace operator. 

a) Let us find eigenfunctions ip(x,y) of the transformed Hamiltonian Hi. On 
M 2 , eigenfunctions of Laplace operator can be chosen in the form of plane waves 
exp(ikr) : 

ip(x, y) = exp(ik x x) exp(ik y y). (9) 

It is clear that any linear combination of eigenfunctions, for which the value k 2 + k 2 
is the same, will also be an eigenfunction. Let us find k x and k y from the boundary 
conditions (E]) , (J6j) : 



k x — + 



27TO 2 



if 1 2nni 
ky ~ b + b 



(10) 



where ni,n 2 G Z. 

Eigenvalues of Hi are found from the equation Hiip(x,y) = Eip(x,y) and have 
the form: 



E 



2m 



' b 



2imi 



<f2 

a 



2irn 2 \ 



(11) 



Let us prove that we have found the complete basis of the Hi eigenfunctions. 
If the torus is without phase lagging, then we obtain the regular Fourier series. 
With phase lagging, not quite regular Fourier series is obtained, however it is re- 
ducible to the regular one: 



ip{x,y) = exp 



( 



, fiy 



\ a 
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Here, \I/ is a regular periodic function for which the normal Fourier series is written 
as an expansion using the basis in L 2 for which the completeness has been already 

proved. Multiplication by exp —i (— — h ) in terms of space L 2 is unitary 



transformation. 



a 



b) Let us find the relationship between the fluxes of magnetic field $i,$2 and 
phase lagging ipi,(p%. The flux of magnetic field is: $ = (7) HdS = (f) rotAdS or, 

after a transformation using Stokes theorem: $ = (j)A.d\. Then the fluxes of magnetic 

field inside the torus and through the hole of the doughnut are respectively equal 
to: b 

$x = f A y (x,y)dy, (12) 
Jo 



$ 2 = / A x {x,y)dx. (13) 
Jo 

Recalling expression (j?]) and substituting it into (fT2l) and ffTB"]) : 



$i = /M)-/(x,0), (14) 

$ 2 = /(a,j/)-/(0,j/). (15) 

It is possible to find the differences in values of f(x,y) in these points using expres- 
sion (jSJ) and boundary conditions (j2J), ([3]), ©, (ED- As a result we obtain: 

^ = -w = (16) 
e 7r 

where % = 1,2; $o is the quantum of magnetic flux. It is necessary to note that in 
this case the magnetic flux is not quantized. 



12. 3D Delta function 

Coulomb wave function of the ground state has the form 

<p B (p) = 8ira(i\ip c (r = 0)\ip 2 , |^ c ( r = 0)| 2 = , Lp p = {f + /iV)" 1 

7T 

and satisfies Schroedinger equation in the momentum representation. 

The values of hyperfine splitting of the ground level of hydrogen-like atom with 
accuracy up to a 5 found on the basis of quasipotential built from the diagrams of 
the order of a 2 and higher, are found most easily if to assume 

<p c (p)^W 3 m\<Pc(r = 0)\. (*) 



Mathematical Physics: Problems and Solutions 



45 



(Infrared singularities in the elements of the amplitude of scattering that appear in 
the framework of this approach, are normally eliminated by cutting off the value of 
virtual 3-dimensional momentum.) 

Derive the relation (-^) by proving the following statements: 

(a) 5(x) = lim - a x G R 1 ; (b) ^ % = lim , 9 - 9 . 9 , x G R 1 ; 
a^o ix a 2 + x 2 2x 2 a^o (a 2 + x 2 ) 2 

0[ J) ) & 

(c) 5(p) = - — - , p = \p\,p ^ 0; (d) 5(p) = lim — — , p = \p\,p ^ 0. 
2irp z a^o ir z {a z + p z ) z 

Direction : equations (a) - (d). containing generalized functions, of the class D' 
should be proved on the space of the basic functions of the class T>. 

SOLUTION 

The above statement {^) is based on the fact that the value of the square of 
module of the wave function (present in a matrix element) in coordinate space at 
r = has the order0 of a 3 . 

Let us prove the (a) statement. 

Let x G R 1 ,(f(x) G V^R 1 ) is the test function. Then 

oo 

lim-—— — -,ip(x)) = / dxlim- — — — ^<p{x) = 

a->0 TTX 2 + CL 2 J J a->0 71 X 2 + CL 2 



Um a f dMx± = lim a 2 V {x) 



a-»0 7T ,/ X 2 + a 2 a->0 7T + GT 

= lim-2^^ = v.(0) = (5(x),^(x)). 

The first and the last parts are underlined to focus on equality of the functionals 
hence the generalized functions themselves. 
Derivation of the (c) relationship: 



S(p) = S(p x )S(py)S(p g ) = (2tt)- 3 J e lkp dk = 

7T 27T oc oo 1 

(^y 3 J J J e ikpcos6 k 2 dksm6d<pd9 = W 2 J j e ikpx dxk 2 dk 

0-1 

1 



(27T) 



-2 



ikpx 



ikp 
o 



k 2 dk = \ n f (e ikp - e~ ikp ) kdk = , \= f e lkp kdk 
(2n)Hp J V ' (2n) 2 ip J 



-l o 



11 Description of hyperfhie splitting is based on computation of quantum distributions being 
squaresd by the wave function. From the explicit form of the wave function in the momentum 
representation it turns that the terms of the lowest order in the expansion of a matrix element are 
a 5 . 
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(2n) 2 ipidp J 2ixp 2np 2 

oc 

where the known relation p5'(p) = — 5{p) is utilized (the proof is partial integration). 

To check (d), we consider the test function G T>(B?) and perform the 

following transformations: 



o->0 7T 2 (p 2 + a 2 ) 2 ' / J a^O 7T 2 (p 2 + a 2 ) 

= [ dfl lim / P dp = 

J a^O TX 2 J (p 2 + a 2 ) 2 



(oo oo 
i / ysfae,^ i rp^{-p^f) 
2 J (p 2 + a 2 ) 2 dp -2j (p> 2 + a 2 ) 2 dp 


{let <f>(p,e,<p) = $(-p,e,<p)} 

dfi lim / P ®'^ dp = 
a->0 27r 2 y (p 2 + a 2 ) 2 

a p 2 $(p, ©, 02) 

cZQ lim — -27riRes„ = j a — — - — (residue in the pole of the 2nd order) 

a^O 27T 2 Qo 2 + a 2 ) 2 

a^O 7T p-Ha Op (p + ia) 2 

aiZ lim — lim 



a^O 7T P^«a ^ (P + id?) 2 (j) + id) 3 

f a ^ f2m$(m, 0,^) (m) 2 $'(m, 6, a?) 2(m) 2 $(m, 6, if) 1 

/ ai 2 lim < ; — 1 ; — ; r~ > 

J a^o Ti \ (2m) 2 (2m) 2 (2m) 3 J 



ai f$(m) $'(m) $(m) 
dfi hm — <^ — i-^- + 



a^o 7T 1 2m 4 4m 



The truth of the statement (b) follows particularly from the proven (c) and (d). 
It is also possible to obtain the relationship (b) using the representation of the 
<5-f unction from the statement (a). 

It is possible to prove the relations (a) - (d) in another way. 

Using Sokhotsky's formula 

P 'mS(x), 



x + iO x 
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which is correct in D', it is easy to derive eq. (a): 

la ,1/1 
hm = hm 



a^O 7T a 2 + X 2 a-H) 2m \X — Id X + id 

= — ( P- + i7r5(x)-P- + i7r5(x) ] = 5(x). 
2ni \ x x J 

It is possible to transform eq. (b) into the following form: 

a^O IX ( a 2 _|_ x 2^ 

Transforming the right part of the last equation, we obtain by recalling eq. (a): 

, 2 ax 2 ,2a/ a 2 

hm ^ = hm — - 



a ~>° 7T (a 2 + x 2 ) Ti a 2 + x 2 \ a 2 + x 



26(x)-lim f (2) 
n (a 2 + x 2 ) 



For any function tp(x) G T>, the following relations are true: 
2a 3 \ t f 2a 3 </9(x)dx 



lim — li m i 

tt (a 2 + x 2 ) 2 / "->oy vr(a 2 + x 2 ) 2 

R 

2 /" 09(at)dt 

hm 



a^OTTj (1+t 2 ) 2 

R 

Um 2 / >(at)-y>(0)] dt , 2y?(0) / dt 



\2 ' 



a^OTry (l+t 2 r vr J (1 + t 2 ) 

R R 

The second integral in the last expression can be computed with the help of the 
transformation t = tga, dt = (t 2 + 1) da: 

7T 

+00 2 

dt / , 71 



[1+t 2 ) 2 



cos ada = — . 



The first integral can be estimated in the following way. Let us split the integral 
into the two integrals: 

lim ^l*"-fd» = 

a^onj (1+t 2 ) 2 
Um r 2_ M at)-m] At Um f 2_ M at)-m] At 

a^O J 7T (1+t 2 ) 2 J IT (1+t 2 ) 2 

l-A;A] R\[-A;A] 
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The first integral on the right side of the last equation tends to zero at a — > 0. 
For the upper limit of the second integral, it is possible to use the value of 

2 2M , 4M f dt 



lim / T^dt 

a->0 



Ml+t 2 ) 2 " 71 J (1+t 2 ) 2 ' 

M\[-A;A] K\[-yl;A] 

independent of a. At A — > oo, this integral also converges to zero. Thus 

lim — n,<^| = <p(Q). 

a^o n (a 2 + x 2 ) 2 J 

It follows that 

2a 3 

lim o = Six). 

tt (a 2 + z 2 ) 2 

Recalling this result, eq. (j2J) converges to eq. (JTJ), which proves the validity of 
eq. (b). 

Eq. (c) should take place on the functions from P(M 3 ), i.e. from that, the eq. 
should follow: 

J S(p)<f(p)dp = j ^^ip(p)dp, (3) 

R 3 R 3 

where <p(p) G P(IR 3 ). Conversing the right part of this equality, we obtain 

/ = I ^$v(P)P 2 d P dn = 



28(p) 



( \ 



P^O 



dp = 5(p)ip(p)dp = 0(0) = v?(0). 



\p\=p ) 

(pip) — . In the last line, function <p was evenly 

47T 

\p\=P 

extended to negative p. The integral in the left part of eq. ([3]) also equals ip(0) which 
proves the validity of that eq. ([3]) as well as of eq. (c). 
Equation (d) automatically follows from (b) and (c): 

5 s (p) = = — lim * = lim ^ . 

'' KFJ 2lip 2 2TT a^O n ( a 2 + p 2f a^O ^ ( fl 2 + p 2)2 



Let us obtain the final approximate equation (+) for the wave function of the 

Sip) 
2np 2 



ground state in Coulomb field using the formula 5(p) = o _ - , valid for the case of 
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spherical symmetry and also using the above proved statements. In the momentum 
representation this function takes the form: 

fcip) = %naii | (p c (r = 0) ' 



(p 2 + a 2 ix 2 ) 
and, using the identity (d), we obtain 

ton<Pc(P) = lim . . = 8tt \(p c (r = 0)\6{p). 

Due to the fact that afi <C 1, the function approximately equals its limit at olll — > 0. 

From the form of Coulomb wave function of the ground state it follows that 
the main contribution into the splitting of the energy levels is due to the momentum 
from the range that satisfies the condition p 2 ~ a 2 fi 2 . As a result, expansion of the 
integrand by p/m would be equivalent to an expansion of the whole integral by a 
(under condition that the integral converges). 

13. Heat conduction equation (heat source presents) 

The temperature on the ends of thin regular rod is maintained constant and 
equals zero. Lateral face of the rod is heat-insulated. The frame of reference is defined 
so that coordinate axis x is oriented along the rod, its ends have the coordinates 
x = and x — I. 

Thermal diffusivity coefficient of the material of the rod equals a 2 . 
Find spatial and time distribution T(x,t) (0 ^ x ^ /, t ^ 0) of temperature 
along the rod in two cases: 

(a) at time moment t — the temperature of the rod is constant, T(x, 0) = T 
at < x < I; 

(b) in the center of the rod, point source of intensity Q is switched on at time 
moment t — 0, and T(x, 0) = at < x < I. 

SOLUTION 

Case (a). In this case, the decision function is defined as the solution of heat 
conduction equation 

T t = a 2 T xx , (1) 
satisfying supplementary conditions 

T(0,t)=T{l,t) = 0^0), (2) 

T(x,0)=T (0<x</). (3) 
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First, let us find the eigenvalues and eigenfunctions of equation ([I]) using separa- 
tion of variables (Fourier method). 

Representing the decision function in the form 

T(x,t) =X(x)U(t) 

and substituting the expression into (pQ), we obtain 

XU = a 2 UX". 

Hence 

a- _^ 

X a 2 U 

X" + X 2 X = 0, (4) 

U + a 2 X 2 U = 0, (5) 

where A 2 =const. Solving eq. (j3J) and (jSJ) with consideration of boundary condi- 
tions (j5J), we obtain: 

/2 / 7TTIX \ 

- sin y — - — J (0 ^ x ^ I, orthonormal system of functions), (6) 



f/ n = C n exp^-(™)%J (t>0), (7) 
where n = 1,2, and C n - real coefficients. It is possible to represent the general 

oo 

solution in the form: T(x,t) = "S^ X n (x)U n (t) . Thus 

n=l 

m , . ^ /2 /7rnx\ ( /irna\ 2 \ , , 

T(x,t) = ^C„y 7 sin(- r jexpf-(- r J t J . (8) 

n=l ^ ' 

According to OH]), expansion of the distribution function of initial temperature 
into series by eigenfunctions ([6]) takes the form 



~w \ v—^ ^ /2 /Tin: 
T(x,0) = ^C n y y sin (— 



'2 /7rna; 

n=l 

where 
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Recalling initial condition 



2 i To l sin 



2 I . 

-T — — cos nn + 1 
/ nn 



2/Tn 



nn 



2V21 



-T , n = 2k + l, 



vr(2£; + ir u ' " ' (9) 

0, n = 2A; + 2 (fe = 0, 1,2, ...)■ 

Substituting the expression (jHJ) into (151) , we obtain the solution for the case (a): 



7T " 2k + 1 

fc=0 



exp 



vr(2fc + l)a\ z \ ^ f-K{2k + l)x 



I 



I 



Case (b). Method I . In this case, the decision function is defined as the solution of 
heat conductivity equation with the source, i. e. boundary problem with the density 
of heat generation described by Dirac 5-function: 



T t = a 2 T xx + — 5\x- 1 - 



T(x, 0) = 0, T(0, t) = T(l, t) = 0. (10) 



Here, c is heat capacity of unit length of the thin tod, Qt (x - - intensity of 

heat generation per unit length. 

We shall search for the solution in the form of the sum of a stationary one (u) 
with a non- stationary one (v): 

T(x, t) = u(x) + v(x, t). 

Stationary solution satisfies 



u 3 



Q ( I 
Six-- 



a?c 



u>(0) = = 0. 



And non- stationary solution satisfies 

v t = a 2 v xx , v(0,t) = v(l,t) = 0, v(x, 0) = —u{x). 

The common solution of stationary equation can be presented in the following 
form (see [2]) 

/ 



U) = d 



X -~2 



C 



2- 
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where indeed the first constant can be found from differential equation 

Q 



Q l 

w M = z-o{x - - 

a z c 2 



C-i 



2a? c 



and C*2 from boundary conditions 



exactly 
Finally 



u(0) = u(l) = 0, 



2a 2 c2 a 2 c4 



_g_ 

2a 2 c 



x 

2 



a 2 c4' 



Non- stationary part of the problem is solved using separation of variables, as in 
the (a) case, only with the boundary condition v(x,0) = —u(x). As a result, we 
obtain T(x, t) = u(x) + v(x, t). 

Note, that expressions for u(x) and v(x, t) can be presented as Fourier series 
expansions (see details in [2]): 



LO(X) 



2QI ^ (_i)fc fix(2k + l)x 
2^ — —sin 1 



n 2 a 2 c^ (2k + I) 2 



(11) 



v(x, t) 



2QI 

7r 2 a 2 c ^ {2k + l) 2 



E 



The answer will be 



Ti(2k + 

j ] ex P 



T{x,t) 



7r(2k + l)a 
I 



t . (12) 



2QI 

ir 2 a 2 c 



n(2k + l)x 
I 



1 — exp 



n(2k + l)a 
I 



(13) 



Case (b). Method II . When solving (1T0|) . let us expand T(x,t) into Fourier series 

oo 

r^, \ \ -* _ , , /7vnx\ 
T(x,t)=J2Cn(t) )» 

n=l 

where 

Cn(0 = y y 71 OM) sin (~y~) 
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In terms of C n (t), the equation (UUp will take the form (use the expansion of 
delta-function) 

"Kna\ 2 2Q . fwn\ 

t) c » (t) + ir m (-) • 

The general solution of the eq. (|14p is the sum of the general solution of the 
uniform equation and partial solution of the non-uniform equation, i. e. 

-2 



Cn{t) 



(14) 



/7ma\ 2 
C n {t) = A n exp I - y—j—j t 



2Q /irn\ /una 



Recalling the initial condition (110)) . we obtain 



2 f f 7TTIX \ 

C n (0) = -\ T{x, 0) sin dx = 0. 



(15) 



(16) 



Considering (TTB"]) . from ( 1T5|) it follows that 



una 



2Q . fim\ ( 



Substituting (1TT1) into (1T51) . we obtain 
r» Art 2( 5 • f nn \ (™ a ^ " 



1 — exp 



/"iina\ 2 



V / 



(17) 



(18) 



Taking into account ffTB)) . we can write the solution of eq. f ll0|) in the form 



Because 



n=l 



exp 



r 



Tina 



sin (~r)- (19) 



/Tin 

sin vT 



(_l) fc , n = 2Jfe + l, 
0, n = 2£; + 2, 



where k = 0, 1, 2, the solution (lT9"j) coincides with (Il3p . 

In the limit t — > oo, the solution (lT9"j) would tend to the stationary solution u(x) 
determined by (ITT1) . 



14. Heat conduction equation with nonlinear add-on 



Burgers' equation is a fundamental partial differential equa- 
tion from fluid mechanics and other areas of applied math- 
ematics. It bears the name of the Dutch physicist Johannes 
Martinus Burgers (1895-1981). For a given velocity of a 
fluid u and its viscosity coefficient u, the general form of 
Burgers' equation has the following form: Vt + vv x = vv xx . 
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Show that it could be linearized by substitution 

v = -2u^- In/, 
ox 

and reduced to the heat conductivity equation f t = vf X x- 



SOLUTION 

Substitute the expression 



-2i//.// (1) 



into the Burgers equation. The result of such substitution is that the all derivatives 
in the Burgers equation obtain the form: 



f f f f f 

v— + 2v—-, v x = -2v— + 2i>—, 

f f f f 3 

2 J XXX n J XX J X A J X 

V ~T ~f 2 — r' 

Substituting the expressions ([2]) into the Burgers equation, we obtain 

fxt fx ft ( fxxx fxxfx 



(2) 



f P v / p 

The obtained equation can be transformed in the following way: 

— fli] = v _ (f xx 

dx\f) dx\f 

Then, for / we obtain almost the equation of heat conductivity (or diffusion) : 

f(x,t)t = vf{x,t) xx + F(t)f{x,t), 

where F(t) is an arbitrary time function. If F(t) = 0, we really obtain the heat 
conductivity (or diffusion) equation. 

Short reference . Suppose that in a certain region of space all particles are moving 
along straight lines parallel to the X axis. 

Let us designate v = dx/dt - the projection of the medium velocity (being the 
function of the coordinate of the point x and time t) on the X axis. The equation 
of free one-dimensional motion of incompressible fluid is written in the form: 

v t + vv x = (3) 

and, as seen, is non-linear. It has a solution in the form of traveling waves the front 
of which is becoming more steep with time and as a result the wave breaks. There 
are many examples of breaking waves from which perhaps the most visual would be 
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formation of the white caps on the sea surface at strong acceleration of the waves 
by the wind. 

Of course, waves breaking does not always take place. There are some existing 
factors that stop process of steeping wave fronts. 

One of such factors is viscosity. If we add the viscosity term to the equation ([3]) 
then we obtain the Burgers equation 



Here, v is the viscosity factor. Within this model, it is possible to describe the waves 
in which the competition takes place between the two opposite processes, steeping 
wave fronts due to non-linearity and quenching due to viscosity. As a consequence 
of such competition, stationary motion can appear. 

The point of interest of the Burgers equation is the existence of exact solution 
built by Hopf [H] and Cole [TO]. Transformation leading to linearization of the Burg- 
ers equation (recalled in the statement of this problem) is called in literature as 
Cole-Hopf transformation. 



v t + vv x = vv. 



XX • 



■ 
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1) "Linear- Nonlinear" response. Non-linear Burgers equation vt + vv x = vv xx can 
be linearized using Coal-Hopf transformation 

, = -2,|>/. (1) 

Here, v(x,t) is the solution of Burgers equation, f(x,t) - solution of the heat 
conduction equation f t = vf xx . 

For the Burgers equation, the initial condition is given as: 

oo 

v{x, 0)| t=Q = v (x), J v (x)dx < oo. 

— oo 

(a) Using the transformation (JTJ) for Vo(x), find the corresponding function 
fo(x) initial condition for the heat conduction equation. 
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(b) The general solution to the Cauchy problem for the heat conduction equa- 
tion is known: 



\JAnvt 



{x-yf 



dy, x,yeR\ t > 0. (2) 



4z4 

Using transformation ([I]), obtain the solution v(x,t) of the Burgers equation. 
Hints 

1. To answer the question (b) use the result obtained in (a) for the present 
problem. 

2. It is convenient to express the answer to (b) using the function 

y 

t;y) = J v (x')dx' + ^(x - yf . 
o 

(c) For Burgers equation, the initial condition is given: 

v(x, 0) = v (x) = — — - — — + . 1 , 2 - 

l + [x — by l + (x + by 

Use for this case the solution scheme developed above in paragraphs (a) 
and (b), find system's response v(x, t). Follow and analyse evolution of the 
obtained solution in time. Utilize the 'Mathematica' package. 

2) Harmony of a flute. 

In woodwind and brass musical instruments, the source of the sound is 
the oscillating column of air. In a pipe, the standing waves emerge. Such 
vibrations occur at certain eigen frequencies. 

Oscillations of pressure in a pipe of length L are described by the wave equation 

d 2 p p d 2 p 
dx^ = Jd^' 

where p is the overpressure (relative to the atmospheric), p ~ density of air in 
the pipe, (3 - modulus of volume elasticity, x - coordinate along the pipe axis 
(see fig. 1), t is time. 

A specific solution of that equation is the function 

p(x, t) = (A cos kx + B sin kx) cos cut. 

(a) Find the values of A, B, k, u when both ends of the pipe are open, and 
also the condition 

p(x = L/2,t = 0) =p > 

is met. 
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x 



(b) On the basis of the solution obtained, analyse the time evolution of the gas 
pressure p(x, t) in the pipe. For this purpose, you are encouraged to use 
the graphics features of the 'Mathematica' software package. 

For your information: 

1) air density under normal conditions is p = 1, 29 kg/ m 3 , and the modulus 
of volume elasticity is = 1, 01 x 10 5 Pa; 

2) the length of a flute may vary widely, so for illustration, it is possible to 
choose L=0.5 m. 

(c) Illustrate the obtained solution using sound synthesis features of the 'Math- 
ematica' software package. Stipulate an opportunity to hear the fundamental 
tone and some overtones of the pipe of variable length. How would the tone of 
the pipe depend on the following parameters: po, (3,L7 

3) From the history of LHC: LEP. 

At the end of the XXth century, the colliding beams experiments on electron- 
positron accelerator have been held in CERN. Such collider is known as the 
LEP-collider (Large Electron- Positron) . The detectors (see figure) recording 
collisions of particles with anti-particles were placed at the intersections of 
the colliding beams. 




Mounting the ALEPH detector ALEPH, end view 



Fig. 2 Fig. 3 

To study the collision pattern, it is required not only to find out which par- 
ticles are born but also to measure their characteristics with high precision, 
reconstruct the particles' trajectories, find out their momenta and energies. 
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Such measurements are held with the aid 
of various types of detectors that coaxially 
surround the place of the collision of the 
particles. In the area of magnetic field, cur- 
vature of a trajectory (see fig. 4) enables to 
find out the momenta of the products of a 
reaction. 



Fig. 4. Curvature of the particles 
tracks in a magnetic field 



On the figure, the event of the birth 
of a neutral K°- meson (kaon) is 
shown. The length of its trajectory 
is 0.1542208 m, the momentum equals 
1.197206 ■ 10~ 18 kgm/s (or 2.240160 
GeV/c), the speed of meson is 0.976200c, 
where c is the speed of light in vacuum. Us- 
ing these data, find out intrinsic lifetime of 
K°- meson, its total and kinetic energies 
(in GeV). 




Fig. 5 



(b) 

In a magnetic field with induction B = 1.52 T, K°- meson decays into tt + 
and 7T _ - mesons with momenta 5.143114 • 10~ 19 kg m/s and 7.027504 • 10" 19 kg 
m/s, respectively. Analyse maximum possible value of radii of the circles of 
lateral motion (with respect to B) of 7r ± -mesons. Also, find the angle of their 
divergence. The elementary charge eo = 1.6 • 10 -19 Clmb. 



4) Virial of gravitational collapse. In classical mechanics of systems executing 
finite motion, the following relationship takes place: 

( K ) = -\(J2^-^)- (3) 

i 

Here, (K) is the mean (for sufficiently long time interval) kinetic energy of 
the system of point particles defined by radius vectors and exposed to the 
action of the forces Fj. 
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A planet revolves around the Sun. Interac- 
tion between the planet and the Sun obeys 
the law of universal gravitation. The mass 
of the Sun is much larger than the mass 
of the planet so the heliocentric reference 
frame can be considered inertial. 

(a) 

Obtain the relationship between the mean 
kinetic (K) and mean potential (U) ener- 
gies of the planet directly from the Virial 
Theorem (J3]). 

(b) 

A planet revolves around the Sun along 
the circular orbit of the radius R. Show 
that the kinetic K and potential U ener- 
gies of the planet on its circular orbit are 
related in the following way 



K 



1 



-U. 



(4) 



A planet revolves around the Sun along 
the circular orbit of the radius R. If the 
mass of the Sun would instantly diminish 
by 2 times, what will be the trajectory of 
the planet? What relationship would be 
given by the Virial Theorem in this case? 




Fig. 6 to the problem 4) (a) 




Fig. 7 to problems 4) (b) , (c) 



5) Waves on Moebius strip. A Moebius strip is a rectangular block < x < a, 
< y < b, where points with coordinates (0, y) and (a, b — y) are glued 
together (see fig.). 

For sufficiently large ratio a/b, the Moebius strip can be implemented nearly 
without stretching as a surface with an edge in three-dimensional space. 

Let the oscillations of the surface of Moebius strip be described by the wave 
equation for the function u(x,y,t) 

u tt — Am = 0. 

The edge of the Moebius strip is free, and hence Neuman's boundary condition 
is set (see fig. 8) 

u y (x, 0, t) = u y (x, b, t) = 0. 
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b 



1 1 



x 







a 



Fig. 8 



1) State boundary conditions on the gluing line (points with the coordinates 
(0, y) and (a, y), y £ [0, b]) corresponding to the longitudinal vibrations {u is a 
small displacement along the surface). Find eigen harmonic oscillations as the 
solutions of the wave equation with corresponding boundary conditions. 

2) State boundary conditions on the gluing line (points with the coordinates 
(0,|/) and (a,y), y £ [0,6]) corresponding to the transverse vibrations (u is 
a small displacement perpendicular to the surface). Find eigen harmonic os- 
cillations as the solutions of the wave equation with corresponding boundary 
conditions. 

6) Collapse of a bubble. Smooth 2-dimensional surface without self-intersections 
in 3-dimensional space is topologically equivalent to a sphere. At the initial 
time moment, the surface bounds the volume V. Points of the surface are 
moving with normally oriented variable velocities. At each time moment, the 
projection of the velocity on the internal normal equals the Gauss curvature 
(product of the two main curvatures) of the surface. Let the surface remain 
smooth during the process of the motion, self intersections do not occur. At 
certain time moment, the surface collapses into a point. What time will it take 
the surface to collapse into a point? 

7) Random problem. Let £i and £2 be positive random variables on probability 
space {Q, J 7 , P}, and such that for all real p £ [a, b], < a < b, 



E£f <oo. 



E is denoted as the operator of mathematical expectation value: 




Prove that their distribution functions coincide: 



F^(x) = P{£i ^ x} = P{£ 2 ^ x} = F i2 (x) for all x £ M. 
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8) Maximal domain for a matrix. Find maximal domain in which the Cauchy 
problem 

dU 

U(x J t)\ t=0 =T(x),—\ t=0 = N(x), (5) 
for the system of equations 

U u - AU XX = 0, (6) 

with the matrix A — ( ^ t^]^U=( Ul ^\ has a unique solution for any 
x G (0,1). 



13/' U 2 



9) Abel's Analogue. In 1823 Abel has been working on the generalization of the 
Tautochrone Problem (to find a curve along which a heavy particle moving 
without friction would reach its lowest position for the same time independing 
on its initial position). Abel has reached to the equation 



-(p(x), (7) 



o 



where f(x) is the decision function, (f(x) - given function. 
Solution of the equation has the form 



ip(t)dt 
\Jx -t 



In the present problem, it is offered to find a solution of the trigonometric 
analogue of the equation 

X 

/ 7== = ,0<t<x<^, (8) 
J ysin(x — t) * 



where tp(x) 



a/cos x 



10) Problem of energy decomposition. Let u(x,t) G C 2 (R x [0, oo)) be a solution 
of the Cauchy initial value problem for one-dimensional wave equation 



u tt - a 2 u xx = 



in R x (0, oo), with initial conditions 

u(x, 0) = g(x), Ut(x, 0) = h(x) 
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where g(x), h(x) are finite functions. 

+00 

Kinetic energy K{t) = - j u 2 (x,t)dx. 

—00 
+00 

Potential energy P(t) = — J u 2 x {x,t)dx. 

—00 

Prove that 

a) K(t) + P(t)a 2 is constant for any t. 

b) K(t) = P(t)a 2 for rather large t. 

11) Dirac Problem. When deriving so called «Dirac equation* in relativistic quan- 
tum mechanics, Dirac has been driven by an idea of «square-rooting» from 
a second order differential operator. 

Find out in terms of square operator of the first order: 

a) a wave one- dimensional operator; 

6) Laplace operator in R 2 . 

12) Certain process for a wave equation. 

Some process is simulated by a function u(x,t), that satisfies the initial condi- 
tions 

sin 2 ix x, < x < 1, du 



0, x < and x > 1, dt ^' °^ °' 



u(x, 0) = 

It is known that even part of this function u r (x,t) satisfies the wave equation 



u u - a2 u r xx = 



in half plane t > 0. Odd part of this function u n satisfies the wave equation 

< - b 2 < x = 0. 

Find the distance between x-coordinates at which u(x,T) has minimal values 
at sufficiently large T. 

13) Maximal domain and a square. For the equation 

Uxx + y/yU X y = (9) 

find maximal domain area on the x-y plane, where 

du du 

u(x,x) = f(x), — — = 4>{x), < x < 1. 

ox ay 

Show that this domain can be divided into 3 parts by straight linear cuts from 
which it is possible to make a square block. What will be the area of such 
square block? 



66 



Series «Modern Problems in Mathematical Physics». Special Issue N a 3 



14) Evaluation of the solution of the ultrametric diffusion type of equation with 
fractional derivative. When solving equations of the ultrametric diffusion type 
(such equations are related to describing conformation dynamics of compound 
systems such as biomacro molecules) , the solutions are often represented in the 
form of exponent series. One of such series is presented below: 

oo 

S{t) = ^2a- i E p (-b-HP). 

Here, Er(z) = > — : r is Mittag-Lemer function, < B < 1, and t is 

i^o T (t 3n + !) 

time, S(t) - probability of finding the system in definite state groups, a >1, 
b >1 - certain parameters. 

Study asymptotic behavior of the function S(t) at t — > oo and find its asymp- 
totic evaluation by t-depending elementary functions. 
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